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1.  INTRODUCTION 


Regenerative  liquid  propellant  gun  (RLPG)  technology  is  sufficiently  mature  to  allow  the 
testing  of  the  first  155  mm  liquid  propellant  gun.  In  suppot  of  the  development  of  this  155  mm 
technology  demoastrator,  test  fixtures  in  30  mm  and  105  mm  sizes  have  been  built  arid  fired.  The 
data  from  all  three  fixtures  have  been  extensively  analyzed  to  better  understand  the  regenerative 
liquid  propellant  gun  (RLPG)  process.  In  this  report,  experimental  data  for  the  155  mm  Concept 
VIC  configuration  arc  reviewed,  and  computer  simulations  of  the  155  mm  tests  arc  prc.scntcd.  In 
previous  reports’*  the  30  mm  and  105  mm  gun  fixtures  were  discus.scd. 

The  structure  for  the  modeling  effort  discu.sscd  in  this  report  was  dictated  by  the  criteria  for 
the  transition  of  the  liquid  propellant  (LP)  program  from  the  Ballistic  Research  Laboratory  (BRL) 
to  the  Army  Research,  Development  and  Engineering  Center  (ARDEC).  The  transition  critciia 
specity  that  it  should  be  po.ssiblc  to  demonstrate  agreement  between  model  and  test  data  for 
M)  mm,  105  mm,  and  155  mm  by  matching  mean,  filtered  pressure-time  data  to  within  5%  for 
damper,  rc.scrv'oir,  chamber,  and  bore  pressures  and  by  matching  muzzle  velocity  to  within  2%.^ 

The  modeling  effort  was  broadened  by  conducting  blind  simulations  of  experimental  test  firings 
using  a  calibration  of  model  parameters  based  on  data  from  similar  test  firings. 

This  re[H)rt  documents  tlic  modeling  of  the  155  mm  Concept  VIC  lest  firings.  First,  a 
description  of  the  Concept  VIC  fixture  is  given.  The  interior  ballistic  model  utilized  has  been 
described  in  previous  publications."’  The  modifications  pertinent  to  the  Concept  VIC  design  were 
described  in  detail  in  the  first  paper  of  this  .scries  which  documented  simulations  of  the  30  mm 
Concept  VIC  gun.‘  The  choice  of  input  parameters  for  the  155  mm  is  described.  These  parameters 
are  primarily  determined  from  the  physical  charactcri.stics  of  the  gun.  However,  .some  parameters 
cannot  be  determined  directly  and  arc  chosen  ba.scd  on  empirical  data.  To  meet  the  transition 
criteria,  test  firings  of  2  liter  and  5  liter  propellant  charges  arc  considered.  The  7  liter  test  data 
were  available  after  the  transition  of  the  program.  Thus,  as  a  extension  of  the  modeling  effort,  the 
baseline  parameters  from  the  5  liter  scries  arc  u.scd  to  predict  the  performance  of  a  7  liter  charge. 
The  igniter  used  in  the  7  liter  tests  differs  from  that  in  the  5  liter  test,  and  differences  due  to  the 
ch.mge  in  igniter  are  discus.scd.  In  general,  simulations  of  the  7  liter  tests  are  reasonable. 

The  key  feature,  and  major  .'^ourcc  of  unccnainty,  in  the  regenerative  interior  ballistic  process 
ts  the  injection,  breakup,  and  combustion  o(  the  liquid  propellant,  A  simple  droplet  breakup  and 
combustion  model  has  been  utilized.  Droplet  size  distribution,  and,  thus,  propellant  combustion. 
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show  mailced  sensitivity  to  other  gun  parameters.  This  sensitivity  and  the  related  variation  of 
droplet  size  distributions  are  discussed. 

2.  THE  CONCEPT  VIC  LIQUID  PROPELLANT  GUN 

A  diagram  of  a  generic  VIC  liquid  propellant  gun  is  shown  in  Figure  1.  The  monopropellant 
in  the  liquid  reservoir  is  prepressurized  and  is  located  between  the  control  rod  (inner  piston)  and 
the  injection  piston  (outer  piston).  A  primer  is  ignited  and  injects  hot  gas  into  the  combustion 
chamber.  As  the  chamber  is  prcscurizcd,  the  control  rod  is  pushed  to  the  left,  opening  the  injection 
orifice.  The  outer  piston  follows  the  control  rod.  injecting  the  propellant  from  the  reservoir  into  the 
combustion  chamber. 

TThe  motion  of  the  control  rod  depends  on  the  damper  (also  referred  to  as  the  buffer)  assembly. 
After  the  initial  seal  between  the  injection  and  control  pistons  is  broken,  the  liquid  pressure  has 
very  litUe  effect  on  the  control  rod.  As  the  control  rod  moves  to  the  left,  damper  fluid  is  forced 
between  the  rod  and  the  extension  on  the  transducer  block.  The  control  rod  >  >  composed  of  turned 
sections  which  define  the  damper  vent  area.  Near  the  end  of  stroke,  the  damj  .  •  vent  area  becomes 
small,  limiting  the  flow  of  fluid  from  the  damper  and  bringing  the  control  rod  to  a  gradual  halt. 

The  injection  piston  will  track  the  control  rod  in  response  to  the  differential  pressure  between 
the  liquid  reservoir  and  the  combustion  chamber.  As  the  injection  piston  moves  closer  to  the 
control  rod,  the  vent  area  will  decrease,  the  liquid  pressure  will  increase,  and  the  injection  piston 
will  slow  down.  As  the  piston  moves  further  from  the  rod,  the  vent  area  will  increase,  the  liquid 
pressure  will  decrca.se,  and  the  piston  will  accelerate. 

Tf;,  fixture  utilized  in  the  firings  discu.sscd  in  this  report  was  developed  under  contract*  by  the 
General  Electric  Company,  Tactical  Systems  Department  (GE).  Figures  1-3  show  representative 
drawings  of  the  piston  and  damper  positions  for  the  GE  155  mm  fixture  at  the  initial,  midstroke, 
and  end-of-stroke  positions  of  the  pistons. 

3.  BASIC  ASSUMPTIONS 

The  basic  physical  assumptions  of  the  model  and  the  changes  in  the  governing  equations 
necessary  to  model  a  VIC  gun  fixture  were  discus.scd  in  the  first  report  in  this  series  which  focuses 
on  the  30  mm  fixture.'  The  second  report  of  the  series  addresses  the  105  mm  fixture.^ 
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Figure  1.  A  Concept  VIC  Regenerative  Liquid  Propellant  Gun,  Initial  Position. 


LP  Reservoir 


Figure  2.  A  Concept  VIC  Regenerative  Liquid  Propellant  Gun,  Middle  of  Stroke. 


'oncept  VIC  Regenerative  Liquid  Propellant  Gun,  End  of  Stroke. 


4.  155  MM  CONCEPT  VIC  DATA 


The  155  mm  Concept  VIC  RLPG  was  tested  using  LGPI846  to  evaluate  Gun  I  ballistic 
performance  with  different  propellant  charge  lengths.  Table  1  shows  the  155  mm  shot  numbers 
together  with  a  brief  description  of  the  test  series.  Muzzle  velocity  and  standard  deviation  in 
muzzle  velocity  are  given  as  reported  by  GE.’  Shot  24  is  excluded  from  the  2  liter  repeatability 
series  since  the  liquid  propellant  metering  system  indicated  that  the  initial  charge  contained  66  cubic 
centimeters  less  liquid  propellant  than  other  firings  in  the  series.  However,  the  liquid  reservoir  was 
full,  the  propellant  charge  being  limited  by  piston  stroke.  It  was  hypothesized  that  the  additional 
volume  in  the  reservoir  (66  cc)  was  occupied  by  ullage  or  water  from  the  purge  cycle.  (Hie  fill 
system  is  routinely  flushed  with  water  between  firings  as  a  safety  measure  during  testing.) 

Besides  the  differences  in  propellant  charge  sizes  in  Table  1,  the  damper  region  was  modified 
throughout  the  firings  The  damper  fluid  is  Brayco  783  in  all  cases,  and  the  damper  does  not 
contain  the  1/32  inch  hole  associated  with  the  30  mm  gun.*  Lubricants,  30-weight  and  20-weight 
oils,  were  pumped  during  the  propellant  fill  process  between  the  outer  piston  and  chamber  wall, 
outer  piston  and  transducer  block,  and  inner  piston  and  transducer  block  to  lubricate  the  pistons, 
provide  hydraulic  support  for  the  outer  piston  shaft,  and  provide  a  barrier  to  prevent  propellant  flow 
between  piston  shafts  and  the  transducer  block.  Two  igniters  were  utilized  in  the  test  program: 

(1)  an  external  solid  propellant  igniter  with  an  initiator  similar  to  that  used  in  the  30  mm  and 
105  mm  firings,  and  (2)  a  hybrid  igniter  (also  called  a  two-stage  igniter)  consisting  of  a  small  solid 
propellant  igniter  which  initiates  combustion  of  a  liquid  propellant  "puddle"  placed  directly  in  the 
combustion  chamber.  The  "puddle"  combusts  and  raises  the  chamber  pressure  before  injection  of 
liquid  propellant  from  the  liquid  reservoir  and  is  referred  to  as  a  pre-positioned  charge.  Projectile 
travel  is  approximately  591.8  cm  with  a  tube  length  of  609.04  cm. 

The  experimental  firings  which  are  simulated  in  this  investigation  arc  shot  17  (2  liter 
characterization  test),  shots  27  and  28  (from  the  2  liter  repeatability  series),  shots  48  and  51  (from 
the  5  liter  repeatability  scries),  shot  58  (5  liter  with  hybrid  igniter),  and  shot  65  (7  liter  with  hybrid 
igniter),  shot  65  is  from  a  7  liter  characterization  series  which  was  fired  in  preparation  for  a 
7  liter  repeatability  scries  to  be  performed  at  a  later  date.  All  data  displayed  in  this  report  have 
been  filtered.  Units  on  graphs  arc  centimeters  (cm),  grams  (g),  milliseconds  (ms)  and  megapascals 
(MPa)  unless  noted  otherwise. 
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Table  1.  155  mm  Concept  VIC  Shot  Numbers 


Shot  numbers 

Description 

1-19 

2-liter  characterization,  solid  propellant  igniter 

20-30 

2-llter  repeatability,  solid  propellant  igniter 

Mean  velocity  at  200  inches  of  travel=393.2  m/s 

with  SD=0.44%  (shot  24  excluded) 

31-46 

5-litcr  characterization,  solid  propellant  igniter 

47-56 

5-liter  repeatability,  solid  propellant  igniter 

Mean  velocity  at  200  inches  of  travel=586.2  m/s 

with  SD=0.25% 

57-61 

5-liter  charge,  hybrid  igniter  characterization 

62-65 

7-liter  characterization,  hybrid  igniter 

Table  2  displays  the  primary  test  parameters  as  reported  by  GE.’  The  projectile  is  initially 
offset  approximately  10.16  cm  from  tube  origin  with  approximately  an  additional  7.08  cm  to  the 
obturating  band.  The  igniter  vents  perpendicular  to  the  tube  axis,  just  behind  the  initial  position  of 
the  projectile.  The  initial  chamber  pressure  rise  due  to  the  igniter  alone  is  determined  by  water 
shots,  (i.e.,  test  firings  in  which  the  reservoir  is  filled  with  water).  The  solid  propellant  igniter 
charge  is  Hercules  HC25SS  with  an  M1B1A2  primer;  the  liquid  propellant  igniter  charge  is 
LGP1846.  The  initial  chamber  volume  includes  the  tube  volume  which  results  from  projectile 
offset  and  is  determined  by  GE  from  a  computer  aided  design  (CAD)  drawing  of  the  gun.  The 
contours  of  the  inner  and  outer  pistons  were  altered  after  the  2  liter  characterization  series,  resulting 
in  a  change  in  initial  chamber  volume. 

Pressure  gages  used  in  this  analysis  include:  two  gages  in  the  same  axial  plane  near  the 
forward  end  (tube  end)  of  the  combustion  chamber,  referred  to  as  D-plane  gages;  two  gages  in  the 
combustion  chamber,  located  in  the  same  axial  plane  just  behind  the  initial  position  of  the  outer 
piston,  referred  to  as  E-plane  gages;  two  gages  in  the  liquid  reservoir,  located  in  the  transducer 
block;  two  gages  in  the  damper  referred  to  as  0-plane  gages;  and  barrel  gages  B1  at  26.64  cm  of 
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Table  2.  155  mm  Concept  VIC  Test  Parameters 


Charge 

2  liter 

2  liter 
repeat 

5  liter 
repeat 

5  liter 
(hybrid 
igniter) 

7  liter 
(hybrid 
igniter) 

Shots 

17 

20-30 

47-56 

57-61 

62-65 

Propellant  Volume  (cm’) 

2134 

2233.6 

5210 

5210 

7217 

Initial  Chamber 
&  Free  Tube  Volume  (cm’) 

1(K)57 

10323 

10323 

10323 

10323 

Projectile  Weight  (kg) 

43.2 

43.2 

43.2 

43.2 

43.2 

Total  Igniter  Charge: 

Solid  (g) 

455 

450 

450 

280 

280 

Liquid  (g) 

0 

0 

0 

133 

115.2 

Piston  Stroke  (cm) 

4.315 

4.453 

10.376 

10.376 

14.387 

Shot  Start  Pressure  (MPa) 

4.0 

4.0 

10.0 

10.0 

20.0 

Bore  Resistance  Pressure  (MPa) 

0.7 

0.7 

3.0 

3.0 

3.0 

projectile  travel,  B3  at  74.497  cm  of  projectile  travel,  B5  at  157.013  cm  of  projectile  travel,  and  B7 
at  260.279  cm  of  projectile  travel.  The  D-plane  gages  used  in  this  analysis  are  D33  and  D138, 
where  the  number  corresponds  to  the  angular  location  of  the  gages  in  the  D-plane.  The  E-plane 
gages  used  are  E222  and  E342.  Not  all  barrel  gages  were  activated  during  each  firing.  The 
control  and  injection  piston  motions  are  measured  by  an  external  Optrons,  and  projectile  motion  is 
measured  by  a  15  GHz  radar. 


The  projectile  shot-start  pressure  and  resistance  profile  has  been  a  continual  source  of  concern 
in  the  modeling  effort.  A  standard  Ml 07  projectile  was  used  in  all  test  firings.  A  standard 
155  mm  rifled  tube  was  to  have  been  used.  However,  during  manufacturing,  the  angle  of  the 
forcing  cone  was  incorrectly  machined  to  1/2  of  the  specified  angle.  The  effects  of  this  change  in 
the  forcing  cone  on  the  resistive  forces  for  the  M107  were  originally  felt  to  be  minimal.  Thus,  as 
a  first  approximation,  experimental  measurement  by  ARDEC  of  resistive  forces  for  a  M107 
projectile  in  a  standard  155  mm  tube  were  utilized.*’  These  values  were  found  to  be  poor  fits  to 
the  experimental  projectile  motion  in  LP  gun  firings  as  determined  from  radar  data.  Therefore,  an 
alternate  model  of  the  resistive  forces  during  early  projectile  motion  was  introduced.  A  constant 
pressure  associated  with  the  beginning  of  projectile  motion,  loosely  referred  to  in  this  report  as  the 
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projectile  shot  start  pressure,  was  determined  for  each  set  of  experimental  data  analyzed.  This 
resistive  force  is  applied  over  an  initial  projectile  travel  equal  to  1.5  times  the  length  of  the 
engraving  band  (which  is  2.54  cm).  This  is  a  "rule  of  thumb"  used  in  solid  propellant  modeling.'® 
The  resulting  theoretical  projectile  motion  was  then  compared  to  the  early  experimental  projectile 
travel,  and  in  general,  the  mle  of  thumb  appears  to  be  reasonable. 

Beyond  the  shot  start  region,  the  resistance  pressure  is  initially  assumed  to  have  a  constant 
value  equal  to  1%  of  the  maximum  chamber  pressure,  another  "rule  of  thumb"  from  solid  propellant 
gun  modeling  for  a  rifled  gun  barrel.'® 

The  projectile  travel,  derived  from  radar  measurements  for  each  test  series,  indicates  that  the 
projectile  shot  start  pressure  varies  in  the  firings.  Radar  data  are  used  to  establish  the  shot  start 
pressures  presented  in  Table  2,  which  give  good  matches  between  the  model  simulations  and  the 
radar  measurements  for  early  projectile  travel.  An  attempt  was  made  to  choose  a  single  value  for 
the  down  bore  resistive  pressure  for  all  test  series.  However,  the  radar  data  from  the  5  liter  firings 
indicated  a  substantially  different  projectile  resistance  profile  than  found  in  the  2  liter  firings.  In 
discussions,  GE"  indicated  that  the  gun  tube  was  lubricated  after  each  shot  up  to  shot  46,  and  after 
shot  46  the  tube  was  not  lubricated.  Since  shots  48  and  51  were  used  in  the  simulation  of  5  liter 
tests,  it  is  assumed  that  the  change  in  resistance  pressure,  compared  to  the  2  liter  series,  may  be 
related  to  the  change  in  test  procedure  as  well  as  (possibly)  to  differences  in  the  initial  pressure  rise 
rate  in  the  combustion  chamber.  However,  interestingly,  radar  from  the  7  liter  firing  shot  65 
indicates  that  the  resistance  profile  had  changed  again,  although  the  initial  pressure  rise,  due  to  the 
igniter  and  initial  venting  of  the  propellant,  is  similar  to  shot  58  in  the  5  liter  series. 

Although  more  complicated  resistive  profiles  could  be  devised  for  a  particular  shot,  a  general 
formulation  which  would  be  applicable  over  all  cases  of  interest  is  desired.  Thus,  the  shot  start 
pressure  is  applied  over  the  first  3.81  cm  of  travel,  and  a  constant  resistive  pressure  is  applied  to 
the  remainder  of  the  projectile  travel.  The  early  projectile  motion  varies  widely  as  a  function  of 
both  time  and  pressure  for  experimental  data  obtained  in  tests  with  "identical"  projectiles.  The 
tolerance  in  the  diameter  of  the  obturating  band,  6.214  inches  to  6.220  inches,  in  a  tube  with  fewer 
than  100  shots  is  not  expected  to  cause  the  observed  variation  in  projectile  resistive  forces.  The 
cause  of  the  observed  variation  is  unknown.  The  effect  of  the  resistance  profile  on  interior  ballistic 
simulations  is  examined  below. 
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5.  GE  155  MM  GUN  HXTURE  -  2  LITER  SHOT  17 


Test  shot  17  is  a  2  liter,  155  mm,  Concept  VIC  firing  from  the  initial  characterization  series 
The  experimental  muzzle  velocity  determined  from  radar  measurements  is  424  m/s. 

The  chamber  pressure  for  shot  17  at  two  gage  locations  in  the  D-plane  is  shown  in  Figure  4. 
The  chamber  pressures  measured  by  the  two  gages  ate  separated  with  a  difference  at  peak  pressure 
of  about  7%.  Data  from  one  liquid  pressure  gage  and  one  damper  gage  are  available  for  this  shot. 
Inner  and  outer  piston  motion  were  measured  using  an  Optron.  Recoil  has  been  subtracted  from 
the  piston  travel  curves  shown  using  the  free  recoil  model  described  in  the  earlier  105  mm  report.^ 
Unlike  the  earlier  105  mm  case,  the  recoil  has  been  scaled  to  give  the  measured  piston  travel  in  the 
155  nun  case.  The  shape  of  the  combustion  chamber  pressure  profile  in  Figure  4  is  unusual 
compared  to  previous  30  mm  and  105  mm  (as  well  as  subsequent  155  mm)  Concept  VIC  data, 
since  it  is  a  nearly  linear,  increasing  function  up  to  peak  pressure. 

Most  of  the  parameters  for  the  gun  code  are  based  on  physical  measurements  of  the  gun 
fixture  or  measured  propellant  or  damper  fluid  properties  (see  Appendix  A).  However,  some  of 
the  parameters  cannot  be  predicted,  and  the  values  are  selected  based  on  the  analysis  of 
experimental  data.  One  such  parameter  is  the  propellant  discharge  coefficient,  which  is  derived 
from  the  experimental  chamber  and  liquid  pressures  and  the  experimental  piston  travels  using  the 
assumption  of  steady-state  Bernoulli  flow.  An  inverse  codc’^-  is  utilized  to  derive  values  for  the 
liquid  reservoir  and  damper  discharge  coefficients,  shown  in  Figures  5  and  6,  respectively.  The 
value  of  the  reservoir  discharge  coefficient.  Figure  5,  is  greater  than  1.0  for  a  large  portion  of  the 
firing  cycle.  This  result  is  not  physically  meaningful  and  is  probably  due  to  small  errors  in  the 
measured  piston  travels,  which  lead  to  large  errors  in  the  computed  vent  area  and/or  the  pressures. 
The  damper  discharge  coefficient,  shown  in  Figure  6,  is  close  to  1.0,  a  reasonable  value  considering 
the  relatively  large  exit  area.  Thus,  model  values  of  the  two  discharge  coefficients  are  fixed  at  a 
constant  value  of  0.95. 

The  interior  ballistic  model  is  first  applied  in  an  optional  mode  using  experimental  chamber 
pressure  gage  D138  as  a  boundary  condition.  The  results  arc  shown  in  Figures  7-10.  The 
comparison  with  the  experimental  liquid  pressure.  Figure  7,  is  good  until  late  in  the  firing  cycle. 

The  compari.sons  with  piston  travels,  in  Figure  8,  and  damper  pressure,  in  Figure  9,  are  quite  good 
with  some  minor  differences.  The  projectile  shot  start  pressure  was  initially  chosen  as  4.0  MPa, 
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Barrel  Prc.ssurc.  Ga-tic  R1  -  Round  17  (line).  Barrel  Pressure  From  Model  With 
Chamber  Pressure  DL'IX  (dot). 


uith  a  barrel  rc.sislive  pressure  of  0.7  MPa  based  on  an  analysis  of  .shot  27  radar  data.  The 
re.sulting  eompari.son  with  ihc  lube  gage  Bl.  loealed  at  26.642  cenlimclcrs  of  projcelilc  travel,  is 
shown  in  f  igure  10.  7'he  projeelile  passes  the  gage  position  ai  the  correct  lime  in  ihe  simulation. 
However,  there  is  some  difference  in  the  magnitudes  of  the  experimental  and  predicted  pressures. 
The  pressure  in  the  tube  at  gage  Bl  rises  after  the  gage  is  uncovered,  since  the  chamber  pressure  is 
still  rising  at  this  point  in  time,  and  the  projectile  is  moving  relatively  slowly. 

I'o  determine  the  e fleet  of  accumulation  on  the  ballistic  process,  the  code  is  run  assuming 
instantaneous  combustion  of  the  propellant  as  it  enters  the  combustion  chamber.  The  eompari.son 
ol  esperimental  and  simulated  chamber  pressures  is  shown  in  Figure  11.  Piston  travels  are 
presented  in  f’igure  12.  .As  can  be  seen  in  Figure  11,  combustion  is  fairly  erUcieni  in  this  shot, 
and  the  chamber  pressure  is  not  strongly  alfecled  by  accumulation.  These  pressure  traces  lack 
the  slowly  increasing  pressure  in  the  stanup  regime,  indicative  of  accumulation,  which  is  a 
ch.iracterislic  of  other  Concept  VIC  data  analyzed.  The  igniter  vents  approximately  10.16  cm  from 
lube  origin  and  is  not  expected  to  inlluence  the  enhanced  early  combustion.  However,  the  chamber 
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Figure  11.  Chamber  Pressure  -  Round  17  -  Gaec  D138  (lincl.  Model  With  Instantaneous 
Burning  (dot). 


i  igure  12  Piston  Travel  -  Round  17  (line).  Model  With  Instantaneous  Burning  (dot) 
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pressures  are  substantially  lower  than  in  the  30  mm  and  105  mm  data.  The  lack  of  significant 
accumulation  may  simply  be  due  to  the  low  pressures  and  correspondingly  slower  pistons  and 
injection  rate,  which  allows  the  injected  propellant  increased  time  to  heat  and  combust  The  shape 
of  the  pressure  curve  may  also  be  related  to  the  bum  rate  of  the  propellant.  Some  experimental 
evidence'^  indicates  that  the  bum  rate  of  the  propellant  is  substantially  lower  below  100  MPa  than 
it  is  above  100  MPa.  Since  accumulation  is  not  significant,  the  assumption  of  instantaneous 
combustion  produces  approximately  the  correct  muzzle  velocity,  a  predicted  value  of  429  m/s 
compared  to  the  experimental  velocity  of  424  m/s. 

However,  Figure  11  suggests  that  there  is  a  small  amount  accumulation  as  evidenced  by  the 
slower  rise  to  maximum  pressure  in  the  experiment.  The  simulation  can  be  improved  with  a  finite 
rate  combustion  model.  In  the  combustion  model  it  is  assumed  that  the  propellant  jet 
instantaneously  breaks  into  droplets  as  it  enters  the  combustion  chamber.  The  number  and  size  of 
the  droplets  define  a  surface  area  for  propellant  combustion,  and  all  droplets  bum  according  to  a 
pressure-dependent,  linear  surface  regression  law,  which  is  derived  from  experimental  strand  burner 
data.““^  For  simplicity,  it  is  assumed  that  the  droplets  all  have  the  same  diameter  at  any  given 
time  and  that  the  total  liquid  surface  area  is  conserved  (i.e.,  the  Sauter  mean  diameter  is  utilized). 
The  droplet  diameter  is  input  as  a  function  of  chamber  pressure.  The  droplet  profile  in  Table  3  is 
derived  by  adjusting  the  droplet  profile  until  the  model  results  match  the  experimental  chamber 
pressure.  Table  3  indicates  that  the  average  droplet  diameter  is  a  constant  100  pm  from  0.0  to 
20.0  MPa,  decreases  linearly  from  100  pm  to  75  pm  from  20.0  to  40.0  MPa,  and  remains  at  a 
constant  75  pm  during  the  remainder  of  the  process. 


Table  3.  Round  17  Mean  Droplet  Diameter  Profile 


Chamber  pressure. 

Droplet  diameter. 

MPa 

pm 

0.0 

100 

20.0 

100 

40.0 

75 

The  liquid  accumulation  and  the  droplet  mean  diameter  can  al.so  be  derived  from  the  inverse 
code  based  upon  measured  piston  positions  to  determine  the  vent  area  and  amount  of  propellant 
injected  into  the  chamber.  However,  approximations  must  be  made,  and  the  results  are  less 
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accurate  than  the  computed  discharge  coefficients.  Computed  liquid  accumulations  and  the  mean 
droplet  profiles  upon  which  the  simulations  are  based  are  presented  in  Figures  13  and  14, 
respectively.  There  is  some  discrepancy  at  early  times  when  a  small  amount  of  liquid  propellant 
has  been  injected.  The  droplet  profile  in  Table  3  is  chosen  to  be  as  simple  as  possible  to 
approximate  the  effect  of  accumulation  and,  arbitrarily,  to  be  a  decreasing  function  of  pressure. 

The  mean  diameter  computed  from  experiment  suggests  a  droplet  profile  which  begins  as  a  fine 
spray,  increases  to  a  maximum  droplet  size  and  then  decreases.  Although  physically  more 
appealing  as  die  vent  area  increases  during  the  startup,  litde  liquid  is  involved  in  the  startup,  and 
the  simple  droplet  profile  is  adequate  to  simulate  accumulation. 

A  comparison  of  the  model  using  the  droplet  profile  in  Table  3  with  experimental  data  is 
shown  in  Figures  15-22.  The  model  chamber  pressure  in  Figure  15  is  not  as  linear  as  the 
experiment,  but  remains  within  the  experimental  gage  measurements.  The  model  liquid  pressure  in 
Figure  16  compares  well  with  the  experiment,  but  again  is  not  as  linear.  The  two  piston  travels  in 
Figure  17  show  overall  good  agreement,  with  a  discrepancy  between  the  model  and  experiment  in 
the  early  injecdon  (outer)  piston  motion.  It  is  noted  that,  since  1.0  millimeter  of  injection  piston 
travel  results  in  injection  of  approximately  75  grams  of  liquid  propellant,  a  small  error  can 
significantly  alter  the  computed  liquid  accumulation.  The  damper  pressures  in  Figure  18  are  in 
reasonable  agreement.  In  order  to  match  the  experimental  muzzle  velocity  of  424  m/s,  the  resistive 
pressure  is  lowered  to  0.7  MPa,  with  the  resulting  comparisons  with  the  barrel  gages  shown  in 
Figures  19-22.  The  predicted  timing  of  the  gage  responses  to  the  passage  of  the  projectile  appears 
accurate.  There  is  an  apparent  loss  of  gage  signal  in  Figures  20-22.  There  is  a  discrepancy  in  the 
magnitude  of  the  pressure  (recorded  by  the  gage)  and  the  model  prediction  in  Figure  19.  A  similar 
discrepancy  is  noted  and  discussed  in  the  next  section  on  the  2  liter  repeatability  series. 
Interestingly,  the  projectile  shows  almost  zero  shot  start  pressure.  Discussions  with  GE"  revealed 
that  the  gun  tube  was  lubricated  between  each  firing  in  the  characterization  scries,  which  may 
contribute  to  the  low  shot  start  pressure. 

6.  GE  155  MM  GUN  FIXTURE  -  2  LITER  SHOTS  27,  28 

The  ten  shots  in  the  repeatability  .series  for  the  2  liter  charge  consist  of  .shots  20-23  and  25-30. 
Two  shots,  27  and  18,  arc  compared  to  assess  repeatability.  In  addition,  the  model  is  used  to 
analyze  the  data.  The  mean  velocity  for  the  2  liter  repeatability  series  at  200  inches  of  travel  is 
393.2  m/s,  with  a  standard  deviation  of  0.44%,  excluding  shot  24,  as  discussed  in  Section  4.’  The 
initial  liquid  reservoir  volume  contains  approximately  50  cubic  centimeters  of  oil  injected  during  the 
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propellant  fill  cycle  to  seal  the  space  between  the  transducer  block  and  the  outer  piston.  The  effect 
of  the  oil  on  the  propellant  combustion  and  thermochemistry  is  ignored,  since  GE"  reports  that  the 
oil  is  most  likely  deposited  in  the  rear  of  the  liquid  reservoir  and  does  not  mix  easily  with  the 
propellant.  However,  in  the  model,  we  assume  the  liquid  reservoir  contains  only  propellant  with 
the  volume  chosen  to  be  consistent  with  recorded  piston  travel. 

The  piston  and  control  rod  in  the  damper  were  modified  from  the  shot  17  configuration.  This 
modification  results  in  slightly  different  initial  liquid  volume,  initial  chamber  volume  (as  derived 
from  a  CAD  analysis),  and  piston  travels.  The  primer  mass  is  slightly  less  than  that  in  shot  17, 
but  otherwise  the  primer  is  the  same.  The  values  of  the  key  gun  parameters  are  listed  in  Table  2. 

A  standard  M107  projectile  was  used  in  this  test  series.  The  projectile-tube  interface  was  the  same 
as  in  shot  17. 

A  comparison  of  the  experimental  chamber  and  liquid  pressures  from  shots  27  and  28  are 
shown  in  Figures  23  and  24.  The  gages  display  the  discrepancies  associated  with  much  of  the 
Concept  VIC  data.  The  four  chamber  pressure  gages  in  Figure  23  record  pressures  consistent  in  the 
startup  region,  but  the  pressures  then  diverge  by  about  7%  at  maximum  pressure.  In  shot  28,  both 
a  higher  pressure  and  a  lower  pressure  were  recorded,  bounding  the  shot  27  data  which  are  more 
consistent.  The  four  recorded  liquid  pressure  histories  in  Figure  24  diverge  in  the  startup  region 
and  show  a  difference  of  approximately  9%  at  maximum  pressure.  There  appears  to  be  little 
correlation  between  the  shots  (i.e.,  the  highest  chamber  pressure  is  recorded  for  shot  28,  but  the 
highest  liquid  pressure  is  recorded  for  shot  27).  This  may  indicate  calibration  errors  in  the  chamber 
gages.  Comparisons  of  piston  travel,  damper  pressures,  and  barrel  pressures  at  gage  B1  for  shots 
27  and  28  are  not  presented;  however,  agreement  is  good.  It  is  noted  that  recoil  has  already  been 
subtracted  from  the  piston  travel  curves.  The  results  of  the  free  recoil  model  are  not  as  accurate  as 
the  105  mm  results;^  therefore,  the  recoil  model  was  adjusted  to  give  the  correct  maximum  piston 
travel.’ 

In  an  effort  to  obtain  a  more  accurate  history  of  the  combustion  chamber  pressure,  GE’ 
suggested  that  the  pressure  recorded  at  the  D-p!ane  gage  location  be  supplemented  by  the  pressure 
recorded  by  the  E-plane  combustion  chamber  gages.  The  D-plane  gages  are  located  at  the  front  of 
the  combustion  chamber,  just  before  the  taper  into  the  gun  tube.  Differences  between  these  gages 
and  oner  chamber  gages  have  been  observed  in  much  of  the  other  155  mm  data.  The  differences 
have  been  attributed  to  gage  drift  caused  by  heating  due  to  increased  gas  flow  and  turbulence 
since  the  gage  is  located  in  the  chambrage  of  the  ehamber;  also,  the  gage  holes  were  milled  too 
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large.  The  two  E-plane  gages  are  initially  covered  by  the  outer  piston  and  do  not  provide  a  record 
of  the  early  chamber  history.  However,  they  are  exposed  to  the  combustion  chamber  environment 
soon  after  the  start  of  piston  travel.  In  Figure  25,  the  recorded  pressures  at  the  D  and  E-plane 
gages  for  Shot  28  arc  compared.  A  close  inspection  shows  that  the  recorded  pressures  at  D138  and 
E342  are  high  and  in  relative  agreement,  and  at  D33  and  E222  are  low  and  in  relative  agreement. 
The  spread  between  the  recorded  pressures  is  on  the  order  of  14%.  Since  there  is  no  apparent 
evidence  of  systematic  differences  between  the  D  and  E-plane  gages,  no  attempt  was  made  to 
combine  recorded  pressures.  The  D-planc  gages  are  used  in  subsequent  analyses  to  be  consistent 
with  earlier  simulations. 

The  inverse  code  is  applied  to  the  data  to  determine  values  of  the  discharge  coefficients  for 
the  liquid  reservoir  and  the  damper.  The  derived  discharge  coefficients  from  the  two  shots  using 
gages  D138  and  LP270  arc  shown  in  Figures  26  and  27  for  the  reservoir  and  damper,  respectively. 
Values  above  1.0  arc  not  physically  meaningful  and  indicate  error  in  pressure  mcasuremcn:  cr 
piston  travel.  It  has  been  noted  that  the  pressure  gages  arc  not  in  good  agreement  and  that  the 
measured  piston  travel  was  adjusted  by  the  recoil  model.  Thus,  the  data  in  Figure  26  are  taken  to 
suggest  litdc  flow  loss,  and  the  discharge  coefficient  for  the  liquid  reservoir  is  chosen  to  be  a 
constant  0.95.  Similarly,  Figure  27  suggests  a  discharge  coefficient  of  0.95  for  the  damper. 

Identical  values  were  used  in  Shot  17  simulations. 

The  projectile  resistive  pressure  is  determined  from  radar  data  using  the  chamber  pressure, 
D138  gage,  as  a  boundary  condition  in  Shot  28.  A  shot  start  pressure  of  4.0  MPa  applied  over 
1.5  times  the  length  of  the  engraving  band  with  the  resistive  pressure  of  0.7  MPa  along  the 
remainder  of  the  bore  provides  good  agreement  between  experimental  and  predicted  projectile  travel 
(sec  Figure  28).  The  radar  signal  is  temporarily  lost  and  regained  during  firing,  producing  the 
break  in  the  experimental  data  in  Figure  28.  A  comparison  of  experimental  and  predicted  pressures 
at  the  B1  gun  tube  gage,  the  only  barrel  gage  recorded  for  this  shot,  is  shown  in  Figure  29. 

Initially,  the  model  is  applied  with  the  Shot  27  chamber  pressure,  D138  gage,  as  a  boundary 
condition.  The  results  are  shown  in  Figures  30  and  31.  The  liquid  pressure  is  reproduced 
accurately  in  Figure  30,  but  the  piston  travels  in  Figure  31  arc  in  poor  agreement.  Although  not 
shown,  the  damper  pressure  agreement  is  adequate,  and  the  early  projectile  travel,  as  shown  in 
Figure  28,  is  good.  To  investigate  the  inconsistency  between  the  predicted  and  experimental  pi.ston 
travels  in  shot  27,  tlic  model  is  applied  to  Shot  28  with  chamber  pressure,  D138  gage,  as  a 
boundary  condition.  These  results  are  shown  in  Figures  32  and  33.  The  liquid  pressure  agreement 
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Figure  27.  Rcsult.s  From  Inverse  Code  -  Damper.  Discharge  Coefficient  (line)  -  Round  27 
Discharge  Coefficient  (dot)  -  Round  28.  Based  on  the  Damper  Gage. 
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in  Figure  32  is  go(xl,  and  the  piston  travel  agreement  in  Figure  33  is  mucli  improved.  The  results 
again  indicate  the  uncertainty  in  the  chamix'r  pressure  gage  measurements. 

Instantaneous  burning  is  a.ssumcd  in  tlie  comparisons  of  liquid  and  chamber  pressures  and 
piston  travels,  shown  in  Figures  34  and  35.  respectively,  for  Shot  28.  It  can  be  seen  that 
accumulation  is  an  important  factor  in  the  ballistic  cycle.  The  instantaneous  combustion  assumption 
results  in  a  predicted  maximum  chamber  pressure  significantly  higher  than  that  recorded.  Interestingly, 
the  opposite  situation  is  seen  in  30  mm  and  105  mm  daia'^  where  the  rapid  bum-off  of  the 
accumulated  propellant  together  with  the  combu.stion  of  entering  propellant  prtxluccs  a  higher 
experimental  maximum  chamber  pre.ssure  than  Uie  model  assuming  instantaneous  combustion. 

Although  there  is  less  accumulation  in  Shot  17  titan  in  Shot  28,  the  energy  relea.se  is  spread  over  a 
longer  [x.'riod  of  time  in  Shot  28.  The  maximum  chamber  pressure  is  reached  in  Shot  17 
(Tigure  II)  at  approximately  13  ms  compared  to  approximately  18  ms  in  Shot  28  (F'igure  34).  Over 
tills  longer  time  period,  the  gas  volume  is  increasing  as  a  result  of  projectile  motion.  Thus,  the  slower 


event  and  larger  gas  volume  result  in  a  chamber  pressure  which  is  lower  than  that  obtained 
assuming  instantaneous  combustion.  The  higher  predicted  liquid  pressure  (Figure  34)  and  faster 
predicted  piston  motion  (Figure  35)  simply  reflect  the  higher  predicted  combustion  chamber 
pressure.  The  experimental  data  are  similar  to  those  of  shot  17  in  the  lack  of  the  relatively  long, 
flat  pressure  during  the  startup.  However,  a  comparison  of  chamber  pressures  for  shot  17 
(Figure  11)  and  for  shot  28  (Figure  34)  indicates  significantly  more  accumulation  in  shot  28.  The 
predicted  projectile  velocity  at  200  inches  of  travel  is  405  m/s,  3.1%  higher,  due  to  the  higher 
predicted  pressures,  than  the  experimental  value  of  393  m/s. 

Accumulation  is  modeled  by  developing  a  droplet  profile  for  the  liquid  propellant  in  the 
combustion  chamber  to  match  the  experimental  chamber  pressure  recorded  at  the  gage  D138 
location.  The  resulting  droplet  profile  is  shown  in  Table  4. 


Table  4.  Round  28  Mean  Droplet  Diameter  Profile  Derived 
From  Chamber  Pressure  Recorded  at  Gage  D138 


Chamber  pressure. 

Droplet  diameter. 

MPa 

pm 

0.0 

100 

20.0 

100 

40.0 

100 

A  comparison  of  the  model  predictions,  using  the  droplet  profile  in  Table  4,  with  experimental 
data  is  shown  in  Figures  36-41.  A  comparison  of  the  chamber  and  liquid  pressures  is  shown  in 
Figure  36.  The  agreement  is  good,  although  the  model  does  not  produce  the  experimental  dip  in 
pressure  between  12  and  14  ms.  The  piston  travels  in  Figure  37  show  good  overall  agreement  with 
some  slight  discrepancies  between  the  prediction  and  experiment.  The  computed  damper  pressure  in 
Figure  38  is  high  and  only  qualitatively  correct.  The  inaccuracy  may  reflect  an  incorrect  damper 
exit-area  profile  due  to  the  initial  seating  of  the  bushing.  The  projectile  travels  in  Figure  39  show 
good  agreement  in  the  early  motion.  The  model  accurately  predicts  the  time  of  the  uncovering  of 
gages  B1  and  B3;  however,  the  magnitudes  of  the  pressures  at  these  locations  are  incorrect  (see 
Figures  40  and  41).  The  predicted  projectile  velocity  at  200  inches  of  travel  is  397  m/s,  which  is 
1.0%  high. 
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Unfiltered  data  were  reviewed  for  Shot  28  in  order  to  investigate  the  discrepancies  in  the 
magnitudes  of  the  predicted  and  experimental  barrel  pressures  seen  in  Figures  40  and  41.  A 
comparison  of  unfiltered  data  with  the  filtered  data  showed  good  agreement.  Thus,  data  reduction 
is  not  considered  to  be  a  cause  of  the  discrepancy.  Comparisons  between  the  two  chamber  pressure 
gages,  D138  and  D33,  and  the  two  barrel  gages,  B1  and  B3,  are  shown  in  Figure  42.  B1  is 
located  at  26.642  centimeters  of  projectile  travel,  and  B3  is  located  at  79.497  centimeters  of 
projectile  travel.  It  is  noted  that  the  recorded  maximum  pressure  at  the  gage  B3  location  is  higher 
than  that  at  gage  Bl.  Although  it  is  physically  possible,  due  to  a  rarefaction  wave,  for  a  gage 
further  down  tube  to  temporarily  record  a  higher  pressure  than  a  gage  closer  to  the  chamber,  no 
physical  mechanism  can  be  identified  which  would  maintain  such  a  pressure  difference.  In 
addition,  since  the  chamber  pressure  is  low,  and  the  projectile  is  moving  slowly,  the  model  predicts 
no  noticeable  rarefaction  wave.  The  predicted  chamber  and  barrel  pressures,  shown  in  Figure  43, 
show  a  small  pressure  drop  from  the  chamber  to  the  first  barrel  gage  and  an  imperceptible  drop 
from  the  first  barrel  gage  to  the  second,  physically  reasonable  events.  Thus,  it  appears  that  there  is 
a  calibration  problem  with  one  or  both  barrel  gages. 

To  probe  this  discrepancy  between  the  barrel  gages  further,  a  simulation  of  Shot  28  using  the 
lower  chamber  pressure  recorded  at  D33  is  developed.  The  recoil  model  is  first  applied  to  the  raw 
Optron  data  in  Figure  44  using  chamber  pressure  measurements  from  D138  and  D33.  Both 
calculations  give  a  maximum  piston  travel  which  is  too  small.  However,  the  lower  pressure  results 
in  less  recoil  and,  hence,  longer  apparent  piston  travel.  Since  the  pressure  would  have  to  be 
significantly  lower  than  recorded  in  the  experiment  to  produce  the  actual  piston  travel  with  the 
recoil  model  used,  it  is  more  likely  that  the  recoil  model  is  inadequate.  However,  as  can  be  seen 
in  Figure  44,  the  piston  motion  is  not  highly  dependent  on  the  exact  recoil  until  end  of  stroke. 

This  does  not  appear  to  be  a  major  factor  in  the  analysis,  and  the  GE  supplied  piston  travel  is 
used. 


In  order  to  produce  a  simulation  of  the  lower  chamber  pressure  measurement  at  D33,  the 
constant  droplet  size  is  increased  from  a  constant  100  pm  to  a  constant  250  pm,  as  shown  in 
Table  5.  To  match  the  experimental  pi.ston  travel  curves  using  the  droplet  profile  in  Table  5,  the 
damper  discharge  coefficient  is  changed  from  0.95  to  1.05.  Although  the  discharge  coefficient 
physically  cannot  be  gicatcr  than  1.0,  the  larger  value  in  the  Bernoulli  equation  is  a  simple  way  to 
uniformly  increase  the  damper  area.  (It  was  noted  to  be  in  question  in  the  discussion  above.)  The 
resulting  comparisons  with  the  experiment  are  shown  in  Figures  45-50.  The  predicted  projectile 
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Figure  42.  Chamber  Pressures  and  Barrel  Pressures  -  Round  28.  Gage  0138.  Gage  D33  - 
Round  28  nineV  Barrel  Gage  B1  (dot).  Barrel  Gage  B3  (dash). 


i  igurc  43.  Chamber  Pressures  and  Barrel  Pressures  -  Model.  Chamber  Prc.s.surc  -  Model  nine). 
Barrel  Gage  D1  (dot). _ Barrel  Gage  B3  (dash). 
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Figure  48  Initial  Projectile  Travel  -  Round  28  (line).  Model  With  Droplet  Burning  - 
Based  on  D138  (dot).  Based  on  D33  (dash). 


Figure  49.  Barrel  Pressure.  Gage  B1  -  Round  28  (line).  Model  With  Droplet  Burning 
Based  on  D138  (dotl.  Based  on  D33  (dashl. 
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Figure  50.  Barrel  Pressure.  Gage  B:^  -  Round  28  (line).  Model  With  Droplet  Burning  - 
Based  on  D138  (doO.  Based  on  D33  (dash). 

velocity  at  200  inches  of  travel  is  now  392  m/s  compared  to  an  experimental  velocity  of  393  m/s. 
The  previous  predicted  velocity  was  397  m/s  using  the  previous  droplet  size  profiles  cited  in 
Tabic  4.  The  chamber  and  liquid  pressures  in  Figure  45  are  in  reasonable  agreement.  Similarly, 
the  piston  travels  in  Figure  46  compare  well.  The  agreement  between  damper  pressures  in 
Figure  47  is  somewhat  improved  compared  to  Figure  38.  However,  the  agreement  of  early 
projectile  motions  in  Figure  48  (using  Table  5)  is  somewhat  poorer  than  obtained  in  the  previous 
simulation  (using  Table  4).  Comparisons  of  experimental  and  predicted  pressures  at  the  B1  and  B3 
locations  using  both  droplet  size  profiles  are  presented  in  Figures  49  and  50,  respectively. 

Figure  49  shows  that  in  the  simulation,  using  Table  5,  the  projectile  arrives  late  at  the  B1  location, 
and  the  resulting  pressure-time  curve  is  qualitatively  different  than  the  experimental  curve  in  shape. 
The  quantitative  discrepancy  is  also  not  improved.  Similarly,  in  the  simulation,  the  projectile 
arrives  late  at  the  B3  location.  It  appears  tliat,  if  the  amplitude  of  the  experimental  barrel  pressure 
were  matched  exactly,  the  predicted  projectile  velocity  would  be  lower  than  experimentally 
observed.  Thus,  it  appears  that  the  pressures  recorded  at  the  barrel  gage  locations  are  in  error,  and 
die  droplet  size  profile  in  Table  4  provides  the  bc.st  overall  simulation  of  the  experimental  data. 
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Table  5.  Round  28  Mean  Droplet  Diameter  Profile  Derived 
From  Chamber  Pressure  Recorded  at  Gage  D33 

Chamber  pressure. 

Droplet  diameter. 

MPa 

pm 

0.0 

250 

20.0 

250 

40.0 

250 

7.  GE  155  MM  GUN  HXTURE  -  5  LITER  SHOTS  48.  51 

The  5  liter  repeatability  series  consists  of  shots  47-56.  Two  of  these  firings,  shots  48  and  51, 
are  analyzed.  The  mean  velocity  reported  by  GE’  at  200  inches  of  travel  for  this  series  is 
586.2  m/s,  with  a  standard  deviation  of  0.25%.  The  liquid  volume  of  5205  cubic  centimeters  used 
in  the  simulation  to  match  the  reported  piston  stroke  (varies  slightly  from  GE  reported  liquid 
volume)  includes  approximately  84  cubic  centimeters  of  oil,  which  is  pumped  during  the  propellant 
fill  cycle  to  seal  the  passage  between  the  transducer  block  and  the  outer  piston.  The  oil  is 
expected  to  remain  segregated  from  the  propellant,  and  thus,  the  properties  of  the  liquid  propellant 
are  not  expected  to  change.  The  oil  volume  is  included  as  propellant.  The  tube  was  not  lubricated 
after  shot  46,  and  it  is  hypothesized  that  increases  in  the  projectile  shot  start  and  resistive  pressures 
shown  in  Table  2  are  a  result  of  the  change  in  the  tube  lubrication  procedure. 

Comparisons  of  the  experimental  data  from  shots  48  and  51  are  shown  in  Figures  51-56. 
Overall,  the  data  are  quite  reproducible.  The  pi.sion  travels  in  Figure  51  show  excellent  agreement. 
However,  recoil  has  already  been  subtracted  from  this  data.  The  recoil  model  was  calibrated  to 
give  the  correct  maximum  piston  travel,  and  the  actual  comparison  between  piston  travels  in  these 
experiments  is  obtainable  only  by  direct  observation  of  the  raw  data.  The  agreement  among  liquid 
pressures  in  the  two  experiments  (two  gages  in  the  same  plane  in  the  liquid  reservoir)  shown  in 
Figure  52  is  excellent.  The  recorded  chamber  pressures  in  the  D-plane  in  Figure  53  show 
discrepancies,  particularly  at  peak  pressure  and  afterwards.  A  close  observation  of  Figure  53  shows 
that,  prior  to  peak  pressure,  the  recorded  pressure  is  higher  at  one  gage  location  in  each  shot. 
However,  at  peak  pressure,  both  D-planc  gages  in  shot  51  record  pressures  which  are  higher  than 
shot  48  by  approximately  5%.  It  is  not  clear  which  pre.ssure  data  arc  more  accurate. 
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Figure  51.  Piston  Travel  -  Round  48  Hincl.  Round  51  (dotl. 


Figure  52.  Experimental  Liquid  Pressure.  Round  48  -  Gage  LP27()  (line)  and  Gage  l.P9n  (doi). 
Round  51  -  Gages  [.P27n  (dash)  and  Gage  Ll^H)  (dot-dash). 
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Figure  53.  Expcrin 
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Figure  55.  Experimental  Damper  Pressure.  Round  48  -  Gage  090  (line)  and  Gage  0270 
(dot).  Round  51  -  Gages  090  (dash')  and  Gage  0270  (dot-dashV 


Figure  56.  Initial  Proicclilc  Travel  -  Round  48  (line).  Round  51  (dot). 
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Figure  54  compares  the  pressures  recorded  by  the  D  and  E-gages  for  Shot  51.  A  close 
inspection  shows  that  the  D138  and  E342  pressures  are  high  and  in  relative  agreement,  and  that  the 
D33  and  E222  pressures  are  low  and  in  relative  agreement.  The  spread  between  the  two  sets  of 
recorded  pressures  is  on  the  order  of  5-10%.  Again,  since  there  is  no  evidence  of  systematic 
differences  between  the  D  and  E-plane  pressures,  the  D-plane  data  are  used  in  subsequent  analyses, 
if  only  to  be  consistent  with  earlier  analyses. 

The  comparison  of  the  recorded  damper  pressures  in  Shots  48  and  51  shown  in  Figure  55 
shows  excellent  agreement.  The  projectile  travels  (derived  from  radar  data),  as  shown  in  Figure  56, 
are  in  good  agreement.  Overall,  the  ballistic  cycle  appears  reproducible,  with  differences  in  the 
measured  combustion  chamber  pressures  possibly  due  to  gage  inaccuracy. 

As  a  first  approximation,  the  experimental  chamber  pressure  recorded  at  the  D138  location  is 
used  as  a  boundary  condition  with  a  value  of  0.95  for  the  discharge  coefficients  for  the  reservoir 
and  the  damper,  as  in  the  2  liter  experiments.  The  resulting  comparison  of  predicted  and 
experimental  liquid  pressures  is  shown  in  Figure  57  for  Shot  48  and  in  Figure  58  for  Shot  51.  The 
agreement  in  Figure  57  is  poor,  while  in  Figure  58  it  is  good.  Due  to  the  inconsistency  of  5-10% 
in  the  experimental  chamber  pressure  data,  these  results  arc  interpreted  as  indication  that  the 
pressure  measurement  in  Shot  51  at  the  D138  location  is  more  accurate.  Therefore,  the  comparison 
of  the  model  prediction  and  experiment  for  the  5  liter  repeatability  series  is  based  on  Shot  51. 

In  order  to  model  the  effects  of  propellant  accumulation,  the  droplet  profile  shown  in  Table  6 
is  developed,  using  the  pressure  recorded  for  Shot  51  at  the  D138  location.  The  droplet  sizes  are 
not  significantly  different  than  those  in  the  2  liter  repeatability  droplet  profile  in  Table  5  up  to 
40  MPa,  but  are  substantially  larger  than  those  in  Table  4,  which  gave  the  best  overall  agreement 
with  the  experimental  data.  The  higher  chamber  pressure  in  the  5  liter  shots  requires  an  extended 
profile.  Comparisons  of  the  model  predictions  using  Table  6  with  experimental  data  from  Shot  51 
arc  presented  in  Figures  59-63.  The  agreement  of  chamber  and  liquid  pressures  in  Figure  59  is 
good,  while  the  piston  travels  in  Figure  60  .show  some  deviation.  The  agreement  of  damper 
pressures  in  Figure  61  is  reasonable,  except  at  peak  pressure  where  a  difference  of  approximately 
15%  is  .seen.  The  experimental  and  predicted  projectile  travel  curves  in  Figure  62  deviate  at  the 
beginning  of  motion;  but  lire  comparison  in  Figure  63  of  pressure  curves  at  the  B1  gage  location, 
the  only  recorded  barrel  gage  in  this  test,  shows  good  agreement  both  in  time  of  projectile  arrival 
and  magnitude  of  the  pressures.  It  is  again  a.ssumcd  that  the  shot  start  pressure  acts  over  1.5  times 
the  length  of  the  obturating  band.  Overall,  the  comparisons  between  the  model  predictions  and 
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experiment  show  more  discrepancies  than  earlier  efforts.  The  major  difference  is  in  the  predicted 
and  experimental  damper  maximum  pressures. 

In  order  to  evaluate  possible  causes  of  the  differences  between  the  model  and  experimental 
damper  pressures  in  Shot  51,  a  number  of  hypotheses  are  examined.  The  first  involves  the 
calculation  of  the  damper  vent  area.  In  the  model,  the  vent  area  is  calculated  as  the  minimum  flow 
area  between  the  bushing  on  the  transducer  block  and  the  contoured  damper  rod,  as  shown  in 
Figure  64.  Thus,  the  location  of  the  defining  point  on  the  bushing  is  initially  taken  as  point  A. 

As  the  rod  moves  to  the  left  (the  bushing  is  stationary),  the  defining  point  on  the  rod  is  relocated 
to  point  B  due  to  the  curvature  of  the  rod.  The  changeover  location  occurs  near  the  maximum 
vent  area  (minimum  bolt  radius)  when  the  areas  M  points  A  and  B  are  the  same.  At  that  time,  the 
location  of  the  maximum  vent  area  is  assumed  to  instantaneously  move  the  0.16  cm  from 
point  A  to  point  B.  This  assumption  is  of  interest  since  this  transition  occurs  near  peak  damper 
pressure,  and  the  true  maximum  vent  area  which  occurs  between  points  A  and  B  is  not  reached  in 
the  model  with  the  redefinition  of  the  vent.  Removing  the  assumption  and  defining  the  flow  area 
only  in  terms  of  point  A,  however,  produced  only  minor  changes  in  the  damper  pressure.  Thus, 
the  method  of  calculation  of  the  vent  area  does  not  explain  the  difference  between  the  model 
prediction  and  experiment. 

As  a  second  approach,  it  was  hypothesized  that  the  damper  control  rod  geometry  is  inaccurate. 
The  computer  model  is  applied  with  the  recorded  chamber  pressure  at  the  D138  location  in  Shot  51 
as  a  boundary  condition.  The  damper  exit  area  near  peak  damper  pressure  is  then  adjusted  to 
attempt  to  match  the  experimental  data.  Figure  65  is  a  comparison  of  the  vent  area  computed  from 
the  engineering  drawings  and  that  required  to  match  the  experimental  damper  pressure.  Figure  66 
shows  the  effect  on  the  calculated  damper  pressure.  The  results  in  Figures  65  and  66  show  that  the 
vent  area  must  be  increased  substantially,  and  as  a  result,  the  predicted  end  of  stroke  occurs  too 
soon.  If  the  experimental  damper  pressure  is  matched  from  the  beginning  of  the  firing  cycle  by 
adjusting  the  entire  vent  area  profile,  the  result  is  a  predicted  end  of  stroke  which  occurs  much  too 
soon.  If  the  experimental  damper  pressure  is  matched  from  the  beginning  of  the  firing  cycle  by 
adjusting  the  entire  vent  area  profile,  the  result  is  a  predicted  end  of  stroke  which  occurs  much  too 
soon,  as  shown  in  Figure  67.  Thus,  tlie  maximum  predicted  damper  pressure  can  be  lowered  to 
agree  with  the  experiment  by  increasing  the  damper  vent  area,  but  this  change  adversely  impacts 
the  piston  motion.  It  appears  that,  if  the  recorded  chamber  pressure  is  accurate,  the  recorded 
damper  pressure  may  be  too  high.  It  is  al.so  possible  that  the  model  is  not  correctly  capturing  the 
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complicated  feedback  between  the  damper  and  chamber  pressures.  The  chamber  pressure  is  found 
to  be  consistent  with  liquid  pressure  in  Figure  58.  Although  the  predicted  damper  pressure 
typically  deviates  from  the  experiment  more  than  other  measurements,  damper  pressure  predictions 
do  not  exhibit  as  large  a  pressure  deviation  in  any  other  shot  from  the  30  mm,  105  mm,  and 
155  mm  firings. 

As  another  hypothesis,  it  is  assumed  that  the  damper  pressure  gage  is  miscalibrated.  The 
result  of  multiplying  the  experimental  damper  pressure  by  1.15  and  then  comparing  it  to  the  model 
prediction  is  shown  in  Figure  68.  The  agreement  is  now  much  improved.  However,  no  firm 
conclusion  can  be  drawn  from  this  result. 

As  a  final  hypothesis,  GE‘*  speculates  that  the  higher  calculated  peak  damper  pressure  is 
attributable  to  excluding  damper  compliance  in  the  damper  pressure  calculation.  As  the  damper 
wall  expands  at  high  pressures,  the  measured  damper  pressure  lowers.  In  the  GE  code,  an  effective 
bulk  modulus  is  calculated  to  account  for  the  damper  compliance.  Our  analysis  has  indicated  that 
this  correction  is  small  and  is  not  sufficient  to  explain  the  difference  between  experimental  and 
model  damper  pressures. 

In  summary,  no  satisfactory  explanation  was  found  for  the  difference  in  predicted  and 
experimental  damper  pressures  shown  in  Figure  61.  It  has  been  assumed  that  simple  Bernoulli 
flow,  with  a  constant  value  of  0,95  for  the  discharge  coefficient,  sufficient  to  describe  the  flow 
in  the  damper  vent.  This  hypothesis  has  proven  reasonable  for  similar  configurations.  It  is 
possible  that  a  complex  discharge  coefficient  model  is  required  due  to  the  sharp  comer  on  the 
bushing  and  to  the  contour  of  the  rod.  (Note  that  a  variable  discharge  coefficient  only  at  maximum 
pressure  would  require  a  value  greater  than  1.0  to  match  experimental  data,  since  the  model 
pressure  is  too  high.)  Although  the  plateau  in  the  experimental  data  near  peak  damper  pressure 
appears  physically  meaningful,  no  mcchani.sm  is  identified  to  match  the  recorded  experimental 
pressure  curve  with  the  given  geometric  description. 

8.  GE  155  MM  GUN  FIXTURE  -  5  LITER  shot  58 


In  order  to  evaluate  the  effectiveness  of  a  hybrid  igniter,  the  155  mm  Concept  VIC  fixture  was 
employed  in  the  configuration  used  in  the  5  liter  ropcatability  scries.  An  external  solid  propellant 
igniter  was  used  in  all  Concept  VIC  firings  before  .shot  57.  In  preparation  for  a  totally  liquid 
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propellant  igniter,  an  intermediate  hybrid  (also  called  a  two-stage)  igniter  was  tested.  The  hybrid 
igniter  consists  of  a  small  external  solid  propellant  charge,  which  vents  into  the  gun  tube  just 
behind  the  initial  position  of  the  projectile,  and  an  unrestrained  liquid  propellant  puddle  (also  called 
a  pre-positioned  charge)  in  the  bottom  of  the  combustion  chamber.  The  pre-positioned  charge  is 
introduced  into  the  chamber  through  the  fill  lines  in  the  transducer  block  just  before  the  pistons  arc 
scaled  and  the  liquid  reservoir  is  filled.  The  liquid  igniter  fill  has  the  added  benefit  of  flushing  the 
liquid  reservoir.  Other  than  the  change  in  the  igniter,  the  gun  configuration  is  identical  to  that  used 
in  5  liter  repeatability  scries. 

The  igniter  model  implemented  in  tlic  code  is  simple.  An  option  is  added  to  the  gun  code  to 
allow  liquid  droplets  to  be  initially  introduced  into  the  chamber.  The  size  of  the  liquid  droplets 
and  the  burning  rate  of  the  propellant  determine  the  initial  pressure  rise  rate.  The  burning  rate  and 
the  thermochemistry  of  the  "igniter"  droplets  arc  taken  to  be  the  same  as  those  of  the  propellant 
injected  from  the  rc.scrvoir. 
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A  comparison  of  chamber  pressures  for  shot  51  using  the  solid  propellant  igniter  and  shot  58 
using  the  hybrid  igniter  is  shown  in  Figure  69.  The  gage  D138  record  was  lost  at  about  10  ms  in 
shot  58.  However,  the  initial  pressure  rise  in  shot  58  is  much  more  rapid  than  in  shot  51, 
apparently  due  to  the  energy  release  directly  in  the  chamber.  The  totally  solid  propellant  igniter 
vents  through  a  long,  narrow  passage;  and  igniter  gases  experience  high  heat  loss.  The  rate  of 
pressure  rise  with  the  hybrid  igniter  indicates  fast,  efficient  combustion  of  the  liquid  puddle.  The 
input  parameters  for  the  model  for  shot  58  are  identical  to  those  of  shot  51,  except  for  the  solid 
propellant  primer  and  the  droplet  profile.  The  droplet  profile  is  modified  as  a  result  of  propellant 
being  injected  from  the  reservoir  into  a  hotter  and  higher  pressure  chamber  environment.  A 
comparison  between  the  droplet  profiles  is  shown  in  Table  7  and  Figure  70.  Since  the  droplet 
profiles  are  strongly  affected  by  small  changes  in  the  igniter,  damper,  and  projectile  shot  start 
pressure,  the  droplet  sizes  should  not  be  interpreted  literally.  However,  the  comparison  indicates 
more  efficient  combustion  in  shot  58,  reflecting  the  differences  in  chamber  pressures. 


Table  7.  Round  58  and  Round  51  Mean  Droplet  Diameter  Profiles 


Round  51 

Round  58 

Chamber  pressure. 

Droplet  diameter. 

Chamber  pressure. 

Droplet  diameter. 

MPa 

pm 

MPa 

pm 

0.0 

300 

0.0 

300 

30.0 

300 

25.0 

300 

40.0 

200 

50.0 

100 

50.0 

200 

75.0 

75 

60.0 

100 

100.0 

20 

70.0 

100 

125  0 

10 

100.0 

25 

150.0 

1 

Comparisons  of  model  simulation  and  experimental  data  for  shot  58  are  shown  in 
Figures  71-74.  The  chamber  and  liquid  pressures  in  Figure  71  are  in  good  agreement.  The 
predicted  control  piston  travel  in  Figure  72  shows  a  comparatively  large  deviation  from  tlie 
experimental  travel,  although  the  injection  piston  travels  agree  well.  Errors  in  either  the  chamber 
pr'-ssure  or  the  predicted  damper  pressure  are  possible.  The  comparison  of  experimental  and 
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predicted  damper  pressures  is  shown  in  Figure  73.  Not  unexpectedly,  the  predicted  damper  pressure 
is  high  compared  to  the  experiment  at  peak  pressure  as  in  Shot  51.  The  early  predicted  and 
experimental  projectile  travels  shown  in  Figure  74  are  in  reasonable  agreement.  The  experimental 
velocity  at  200  inches  is  596  m/s,  vs.  a  predicted  velocity  of  597  m/s. 

The  hybrid  igniter  appears  to  provide  enhanced  gun  ignition  and  more  efficient  combustion  of 
the  injected  liquid  propellant.  A  simple  igniter  model  is  sufficient  to  model  the  hybrid  igniter. 

The  resulting  comparisons  between  the  model  prediction  and  experiment  are  similar  to  previous 
results. 

9.  GE  155  MM  GUN  HXTURE  -  7  LITER  SHOT  65 

Shot  65  is  part  of  a  7  titer  characterization  series  consisting  of  rounds  62-65.  The  gun 
configuration  is  similar  to  that  used  in  the  5  liter  firings.  The  transducer  block  is  moved  rearward 
to  accommodate  the  additional  liquid  in  the  reservoir,  resulting  in  a  longer  piston  travel.  The 
hybrid  igniter  is  used  with  slightly  less  liquid  propellant  in  the  pre-positioned  charge  (liquid  puddle) 
than  for  Shot  58.  The  initial  chamber  volume  and  pistons  are  unchanged  from  the  5  liter 
repeatability  series.  The  damper  rod  is  extended  by  inserting  a  1.75  inch  straight  section  near  the 
middle  of  the  piston  stroke  to  accommodate  the  longer  piston  travel. 

The  droplet  profile  shown  in  Table  8  is  developed  based  on  the  experimental  chamber  pressure 
recorded  at  the  D33  location,  since  the  signal  from  pressure  gage  D138  was  apparently  lost  later  in 
the  ballistic  cycle.  Initial  analyses,  based  on  the  experimental  chamber  pressure  and  projectile 
di.splacement,  indicated  that  the  shot  start  pressure  must  be  doubled  from  10.0  MPa  in  the  5  liter 
test  series  to  20.0  MPa  in  the  7  liter  scries,  even  though  the  projectiles  are  all  M107  projectiles, 
and  the  initial  chamber  pressure  rise  rates  in  the  two  experiments  are  similar.  As  before,  the  shot 
start  pressure  is  applied  over  1.5  times  the  length  of  the  obturating  band,  a  total  of  3.81  cm  to 
obtain  reasonable  comparison  with  the  experimental  data.  This  will  be  discussed  in  more  detail 
below.  The  comparisons  between  the  model  predictions  and  experiment  are  presented  in 
Figures  75-80. 

The  comparison  of  predicted  and  experimental  chamber  pressures  in  Figure  75  shows  good 
agreement  up  to  peak  pressure,  with  tlie  usual  deviation  during  the  expansion  phase.  The  model 
does  not  predict  the  llattcned  chamber  pressure  curve  near  peak  pressure.  The  gage  with  the  higher 
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Table  8.  Round  65  Mean  Droplet  Diameter  Profile  Derived 
From  Chamber  Pressure  Recorded  at  Gage  D33 


Chamber  pressure, 
MPa 


0.0 

25.0 

50.0 

75.0 

100.0 

125.0 

150.0 


Droplet  diameter, 
pm 


300 

300 

125 

100 

60 

30 

15 


recorded  pressure,  D138,  lost  the  pressure  signal  during  the  ballistic  cycle.  The  agreement  between 
predicted  and  experimental  liquid  pressures  in  Figure  76  is  excellent.  The  injection  piston  travels. 
Figure  77,  are  in  good  agreement;  however,  there  the  control  piston  travels  show  some  deviation. 
The  same  characteristic  was  noted  in  the  5  liter  test  series,  shot  58.  The  predicted  damper  pressure 
in  Figure  78  is  in  good  overall  agreement  with  experiment  with  an  overpressure  at  peak.  Predicted 
and  experimental  projectile  travels.  Figure  79,  are  in  good  agreement,  although  the  early  projectile 
motion  in  the  experiment  occurs  somewhat  sooner.  The  one  recorded  barrel  gage  B3  comparison  in 
Figure  80  is  good.  The  experimental  muzzle  velocity  is  683  m/s  at  200  inches  of  travel,  while  the 
predicted  velocity  is  682  m/s. 

The  flattened  chamber  pressure  curve  in  Figure  75  has  not  been  observed  in  experimental  data 
since  medium  caliber  testing  early  in  the  liquid  propellant  program.  It  is  due  to  the  straight  section 
on  the  damper  rod  which  allows  the  system  to  reach  a  near  steady  state  condition.  The 
regenerative  pistons  pass  the  straight  section  on  the  conuol  rod  just  before  the  plateau  in  the 
chamber  pressure  and  achieve  a  near  steady  state  condition  during  the  travel  over  this  straight 
section.  As  the  damper  vent  area  begins  to  close,  the  piston  momentum  maintains  the  equilibrium 
condition  for  a  short  time.  The  control  piston  then  decelerates  in  response  to  the  decreasing  vent 
area,  and  the  injection  piston  momentum  causes  the  sharp  rise  in  the  liquid  pressure  at 
approximately  15  ms  (see  Figure  76),  As  the  injection  piston  decelerates,  the  injection  area  is 
reduced,  and  chamber  pressure  decreases. 
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The  hybrid  igniter  appeared  to  function  well  in  the  7  liter  test  as  it  did  in  the  5  liter  Shot  58. 
The  only  anomalies  in  the  data  are  the  increased  projectile  shot  start  pressure  and  higher  amplitude 
pressure  oscillations  than  previously  observed,  particularly  in  the  damper.  The  pressure  oscillations 
will  not  be  discussed  in  this  report  and  are  the  subject  of  a  separate,  ongoing  investigation. 

10.  FURTHER  MODELING  -  PREDICTION  OF  7  LITER,  SHOT  65 

Although  the  model  results  presented  above  demonstrate  excellent  agreement  with  the 
experimental  data,  values  of  die  empirical  parameters,  discharge  coefficients  for  the  reservoir,  and 
the  damper  and  the  mean  droplet  diameter  size  distribution  are  derived  from  the  experimental  data. 

A  longer-term  goal  of  the  liquid  propellant  interior  ballistic  modeling  program  is  prediction  of 
performance  for  untested  gun  designs.  Therefore,  it  is  desirable  to  obtain  values  of  these  parameters 
which  are  independent  of  a  particular  fixture.  Thus,  parameters  derived  for  the  5  liter,  Shot  58  are 
used  in  a  simulation  of  the  7  liter,  Shot  65. 

A  constant  value  of  0.95  for  the  discharge  coefficient  for  both  the  liquid  reservoir  and  the 
damper  was  used  in  all  the  155  mm  modeling.  The  droplet  profile  for  the  5  liter  igniter.  Shot  58, 
is  used  since  the  igniter  is  similar  and  vents  into  the  same  volume  in  both  fixtures.  The  prediction 
is  then  not  dependent  on  the  igniter  model.  The  results  arc  shown  in  Figure  81  using  the  droplet 
size  distribution  in  Table  7.  The  initial  chamber  pressures,  which  are  due  primarily  to  the  igniter, 
agree  well  as  expected,  but  the  predicted  and  experimental  pressures  deviate  at  5.0  ms.  The 
subsequent  predicted  pressure  rise  is  too  slow,  and  although  the  maximum  pressures  are  within  5%, 
the  shape  of  the  curve  is  incorrect.  The  predicted  velocity  at  200  inches  of  projectile  travel  is 
660  m/s,  while  the  experimental  velocity  is  683  m/s,  a  difference  of  -3.4%. 

The  inconsistency  can  be  traced  to  the  early  projectile  motion.  Although  the  initial  chamber 
pressure  rise  rate  due  to  the  igniters  is  similar,  as  shown  in  Figure  82,  the  early  projectile 
motion  varies  considerably.  Shown  in  Figures  83  and  84  arc  the  experimental  projectile  travels 
derived  from  radar  for  Shots  58  and  65  vs.  chamber  pressure  and  time,  respectively.  The  two 
projectiles  arc  "identical”  and  yet  behave  quite  differently.  The  model  uses  a  resistive  pressure 
determined  by  applying  a  constant  shot  start  value  over  tlic  first  3.81  cm  of  travel,  the  first  .small 
segment  of  the  curves  in  Figures  83  and  84.  A  constant  shot  start  value  of  10  MPa  simulates  Shot 
58  well,  but  a  shot  .start  pressure  of  20  MPa  is  needed  to  simulate  Shot  65.  It  is  not  apparent  why 
the  resistive  profiles  arc  so  different  between  these  two  shoLs.  The  physical  mechanism  for  the 
di  (Terence  is  not  apparent  and  is  the  subject  of  current  study. 
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However,  if  the  change  in  projectile  shot  start  pressure  is  recognized  by  the  model,  the  result  of 
the  prediction  of  shot  65  from  shot  58  is  much  improved,  as  shown  in  Figure  85.  The  predicted 
and  experimental  chamber  pressures  are  now  in  better  qualitative  agreement.  The  comparison  at 
peak  pressure  is  about  the  same  as  in  Figure  81.  (0.7%).  The  predicted  projectile  velocity  at 
200  inches  of  travel  is  now  687  m/s  compared  to  an  experimental  value  of  683  m/s.  a  difference  of 
+0.6%.  Thus,  the  projectile  resistive  pressure  has  an  impact  on  the  predicted  velocity  and  the  shape 
of  chamber  pressure-time  curve.  While  the  impact  is  not  surprising,  the  deviation  in  projectile 
resistive  pressure  must  be  quantified  in  order  to  achieve  good  predictive  results. 


Figure  81.  Chamber  Pressure  -  Round  65  (linel.  Model  With  Droplet  Burning.  Shot  Start 
Pressure  of  10.0  MPa  (dot). 
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11.  CONCLUSIONS 


A  lumped  parameter  interior  ballistic  model  has  been  developed  for  the  Concept  VIC 
regenerative  liquid  propellant  gun,  incorporating  the  major  features  of  this  configuration.  Most  of 
the  required  parameters  are  physical  dimensions  of  the  gun  and  properties  of  the  propellant,  which 
can  be  chosen  without  reference  to  experimental  firing  data.  However,  the  parameters  which  must 
be  fixed,  based  on  analyses  of  gun  test  data,  are  the  reservoir  and  damper  discharge  coefficients  and 
the  droplet  diameter  profile. 

A  constant  reservoir  discharge  coefficient  of  0.95  appears  reasonable  for  all  the  30  mm, 

105  mm,  and  155  mm  Concept  VIC  firings  studied  thus  far.  The  damper  discharge  coefficient  is 
chosen  to  be  0.95  for  all  155  mm  cases  studied  and  has  varied  from  0.6  to  0.95  in  the  30  mm  and 
105  mm  cases.  It  has  been  possible  to  assume  a  coastant  discharge  coefficient  for  all  the  dampiers 
studied,  although  the  agreement  between  the  model  prediction  and  the  experiment  is  least  satisfactory 
in  the  comparison  of  damper  pressures.  It  is  suggested  that  the  complicated  flow  pattern  may  lead 
to  a  variable  discharge  coefficient.  However,  the  constant  discharge  coefficient  assumption  provides 
good  overall  agreement  between  the  experiment  and  the  model  prediction.  The  mean  droplet 
diameter  also  varies  for  each  gun  configuration.  There  are  no  obvious  correlations  for  the  droplet 
diameter  profiles  for  the  various  gun  configurations,  but  propellant  injection  and  combustion  are 
strongly  influenced  by  the  characteristics  of  the  igniter,  the  damper,  and  the  rcgcnerati'/c  pistons,  as 
well  as  the  chamber  pressure. 

In  general,  tlie  model  can  be  used  to  match  cxperimenial  data  from  155  mm  Concept  VIC 
RLPG  tests,  using  reasonable  values  for  the  empirical  parameters. 

In  summary,  the  following  observations  seem  appropriate: 

(1)  Di.scrcpancies  within  a  given  experimental  data  .set,  particularly  in  chamber  pressure 
measurements  and  even  between  pre.ssurc  gages  in  the  same  plane,  make  accurate 
comparisons  between  the  model  predictions  and  experiment  dirficult.  The  source  of  these 
di.screpancies  is  unknown;  however,  many  gages  show  deviation  early  in  the  ballistic  cycle. 

(2)  Overall,  simulations  of  the  155  mm  Concept  VIC  are  in  good  agreement  with 
experimental  data.  However,  in  order  to  develop  a  predictive  capability,  a  more  general 
combustion  model  is  needed.  Several  approaches  are  currently  under  investigation. 
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(3)  The  damper  submodel  remains  somewhat  unsatisfying  as  in  the  earlier  30  mm  and 

105  mm  modeling  studies.  The  performance  of  the  gun  is  sensitive  to  small  changes  in  the 
damper  since  it  directly  influences  the  motion  of  the  control  piston,  and  therefore,  propellant 
injection  velocity  and  rate  of  propellant  flow  into  the  chamber.  However,  an  improved 
model  may  require  an  empirically  determined,  variable  discharge  coefficient. 

(4)  The  variations  in  the  projectile  resistive  pressure  profile  over  the  first  4  cm  of  motion 
mu.st  be  understood  or  controlled  in  order  to  accurately  predict  gun  performance. 
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APPENDIX  A; 

INPUT  AND  OUTPUT  FILES  -  SHOT  17 
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Below  is  a  listing  of  the  job  stream  for  the  2  liter  characterization  shot  17.  The  numbers  and 
labels  at  the  left  are  read  in  by  the  code.  The  comments  at  the  right  are  for  identification  by  the  user 
and  do  not  affect  the  actual  code.  Following  is  the  summary  shfet  from  the  output  file.  A  description  of 
the  input  and  a  brief  de.scription  of  the  source  for  the  input  variables  are  given  first. 

The  first  line  is  merely  a  label.  It  lists  the  filename  of  the  input  job  stream  and  a  brief  description 
of  the  problem. 

The  initial  offset  of  the  projectile  and  the  total  distance  traveled  by  the  projectile  before  muzzle 
exit  are  given.  The  diameter  of  the  gun  tube  is  given.  The  projectile  and  piston  weights  are  entered 
(measured  by  GE). 

The  initial  volume  of  the  liquid  reservoir  is  given.  GE  reports  in  Table  2  a  volume  of  2134  cm’ 
and  a  piston  stroke  of  1.699  inches  =  4.315  cm.  To  maintain  consistency  in  the  code,  the  liquid  volume 
was  taken  to  be  2166  cm’,  with  a  piston  stroke  of  1.7  inches  =  4.318  cm.  The  initial  free  volume, 
including  the  chamber  and  the  tube  volume  from  the  offset  of  the  projectile,  is  derived  by  GE  from  a 
C.AD  drawing  to  be  10,057  cc.  The  initial  area'  of  the  reservoir  (including  control  rod)  and  chamber  are 
derived  from  the  engineering  drawings. 

The  VENT4  option  is  chosen  (VIC  gun).  The  control  rod  radius  vs.  relative  piston  motion  is 
derived  from  the  drawings.  The  zero  point  is  defined  where  the  outer  piston  and  control  rod  first  fit 
together.  The  positive  direction  is  to  the  left  (direction  of  the  piston  stroke).  In  this  case  only  the  parts 
of  the  control  rod  to  the  right  (negative  direction)  arc  rclevcnt,  since  the  outer  piston  cannot  move  to  the 
left  with  respect  to  the  control  rod.  The  area  of  the  hole  in  the  outer  piston  is  computed  from  the 
engineering  drawings.  There  is  not  a  grease  dyke  on  this  piston. 

The  piston  resistance  is  set  to  zero,  since  this  is  not  usually  large  enough  to  be  important.  The 
discharge  coefficient  into  the  chamber  is  set  at  a  constant  0.95  (sec  discussion  in  report).  The  discharge 
coelficient  into  the  gun  tube  is  .set  equal  to  one  (no  losses). 

The  How  into  the  chamber  is  modeled  as  steady  state  Bernoulli  flow  (FLLJXl).  There  is  only  one 
\eut  hole,  fhe  piston  tliickne.ss  is  irrelevant  for  this  option.  The  How  into  the  gun  tube  is  steady  state 
l'^e^tropk  How  (1  I  t  I  he  shot  start  pressure  is  applied  over  a  distance  etjual  to  the  1.5  times  the 

lengili  of  the  obturating  band,  a  length  o/  1.5  inches  =  .VSl  cm  (see  discussicu)  in  rep<in).  The  value  of 
4  (1  MPa  gives  very  good  agreement  with  the  early  projectile  travel.  After  shot  start,  the  resistance 
pressure  is  tjuickly  lowered  to  0.7  MPa  (see  discussion  in  report).  The  155  mm  is  a  rilled  gun. 
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Next,  the  physical  properties  of  the  propellant  HAN1846  are  given. 


The  liquid  is  pre-pressurized  to  3.45  MPa.  The  initial  chamber  pressure  is  one  atmosphere. 

I  he  droplet  diameter  is  read  in  as  a  function  of  chamber  pressure  (DROP4).  The  burning  rale  is 
read  as  a  two  part  function.  There  is  some  evidence  that  there  is  a  break  in  the  slope  of  the  burning  rate 
just  under  100  MPa.  However,  since  the  rate  at  higher  than  100  MPa  is  not  known,  the  rate  that  was 
measured  for  lower  pressures  is  used.  The  droplet  diameter  table  has  been  chosen  to  match  the  chamber 
pressure. 

The  primer  is  injected  in  the  form  of  hot  gas  (PR1M3).  The  actual  primer  mass  is  455  g.  From 
tltc  water  shots,  the  primer  presurizes  the  chamber  in  about  4.0  ms.  A  heat  loss  factor  of  0.45  is  used. 

I  hal  is,  only  45%  of  the  energy  of  the  primer  actually  pressurizes  the  combustion  chamber.  From  the 
GF  water  shot,  chamber  pressure  reaches  about  18  MPa  in  4.0  ms.  The  above  value  leads  to  18.4  MPa 

at  4  ms. 

The  default  models  for  the  heat  loss  to  the  gun  tube  and  the  air  shock  arc  used.  Heat  loss  to  the 
combustion  chamber  walls  is  ignored. 

The  most  complicated  Lagrange  model  (TUBE4)  is  chosen.  The  model  will  take  into  account  the 
larefaction  wave  after  burnout  of  the  propelltint. 

The  buffer  model  is  chosen  (BUFF2).  The  areas  of  the  buffer  side  of  the  control  rod  and  the  hole 
in  ilie  block  are  from  the  engineering  drawings.  The  discharge  coefficient  is  set  equal  to  0.95  (see 
discussion  in  report).  The  initial  volume  is  estimated  from  the  drawings.  The  buffer  is  originally 
pressuri/ed  to  3.45  MPa  to  reduce  ullage.  The  buffer  fluid  is  Brayco  783,  which  has  a  density  of 
i)  SSS5  gy'ce.  Unfortunately,  the  bulk  modulus  for  the  damper  fluid  has  not  been  measured.  The  bulk 
mi'diilu^  for  a  similar  fluid,  Brayco  756,  has  been  measured  up  to  21  MPa  at  24.4  C.  This  is  fit  by  a 
linear  lunclion  to  obtain  the  derivative  of  the  bulk  modulus  with  rcspeci  to  pressure.  This  means  that  the 
bulk  modulus  of  the  buffer  lluid  is  uncertain,  especially  at  high  pressures. 

The  center  bolt  radius  is  given  as  a  function  of  control  rod  travel.  The  bolt  on  the  back  of  the 
:  ooi  lod  tias  been  turned  to  fonn  this  vent. 

i  lie  code  Will  pnnt  out  results  everv  0  2  ms  (TI.NC).  Because  tlie  code  must  oltcn  change  the  time 
:  ;  .  '  more  clliv.icnt  to  restnet  the  maximum  time  step  (HTOP)  The  error  controls  LP.S  and  SRFf 

.I'e  ;...ven  upual  values. 
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The  integration  method  flag,  MF,  is  set  to  10  (variable  step  size,  variable  order  Adams  method).  For 
the  Concept  VIC  configuration,  this  ran  more  rapidly  than  the  backwards  differentiation  formula 
previously  used.  KWRITE  is  set  to  zero  to  eliminate  diagnostic  messages.  A  time  limit  of  TMAX  is  set. 
If  the  code  takes  longer  than  TMAX  seconds  to  execute,  the  code  will  stop  gracefully  and  write  tiie  usual 
summary  pages. 

The  code  is  only  to  be  integrated  once  (REPl)  and  the  chamber  pressure  will  be  computed  normally 
(CHAMl). 
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dtl7  -  VIC  -  2  liter  Charge  - 
10.16  591.8 

15.5 
43200, 

109000. 

2166.  8140. 

641.19  855.02 

vent4 
11 

-3.0620  3.4925 

-2.2873  7.3762 

-0.0521  7.8105 

0.0000  7.8359 

0.7798  8.2169 

1.0820  8.3185 

1.3360  8.3452 

1.5900  8.3185 

1.7170  8.2868 

3.3388  6.6650 

34.2760  6.6650 

192.898  0.0 

61400 
no  close 


pisl 

2 

0.0  0.0 

1.0  0.0 

disl 

2 

0.0  .95 

1.0  .95 

disl 

2 

0.0  1.0 

487.17  1.0 


Drop  -  CD=0.95  -  CD5=0.95 

offset  proj  travel 

gun  tube  diameter 

proJ  v/eight 

piston  weight 

vl  v3 

al  a3 

moving  central  bolt 

bolt  radius  versus  piston  travel 


ahole  agres 
rodwt 

pistons  initially  meet  at  point 
piston  resistance 


dis.  coeff.  vs.  piston  travel 


dis.  coeff.  vs.  proj  travel  -  tube 
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fluxl 

steady  state  mass  flux  formulation 

1 

1.0 

nvo  pth 

flux2 

isentropic  flow  into  tube 

projl 

proj  resistance 

4 

0.0 

4.0 

3.81 

4.0 

3.82 

0.7 

591.8 

0.7 

1.43 

5350.0  9.11 

rhO  kl  k2 

4035.5 

1.2226 

energy  gamma 

66.9 

.04988 

surface  tension  kinematic  viscosity 

22.848 

.677 

mol  wt  gas  covolumc 

3.45 

0.1  0.0 

pi  p3  epsml 

drop4 

droplets 

0.01 

95.2590 

ddr  pbr 

1.64 

.103 

adrl  bdrl 

1.64 

.103 

adrl  bdrl 

5 

0.0 

.0100 

20.0 

.0100 

40.0 

.0075 

60.0 

.0075 

1000.0 

.0075 

prim3 

inject  hot  gas 

455.0  1 

3.004  0.45 

primer  mass  injection  time  heat 

loss 

0.0 

puddle  mass 

hcall 

no  heat  loss  to  chamber  walls 

hcal2 

heal  loss  to  gun  tube  walls 

300.0 

1.0 

lube  temp  fudge  factor 

shock2 

air  shock 

0.1 

300. 

airp  airl 

1.4 

28.84 

airgam  airmw 

tubc4 

unsteady  lagrangc  distribution 
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bum 

buff2 

139.555  35.04845  0.0 
.95 

898. 

3.45  3.45 

0.8885  1600.  11.3 
0.0  1.4  0.727  0.0 

radius 

6 

0.0000  3.3084 

0.1067  3.3084 

1.0465  3.2576 

1.8288  3.2576 

3.0988  3.2830 

4.3180 
l.OOc-05 
10  0 
20.0 
repl 
chaml 


buffer  -  brayco  783 
apis  ahole  acir 
cd5 
v5 

p5  pout5 

rhO  kl  k2 

cps5  gam5  cv5  b5 

Read  in  damper  rod  radius 

nrat 


3.3269  2.00e-04 
l.OOe-08 


l.OOe-05 
eps  srec 
mf  kwrite 
tmax 

integrate  once 
compute  p3 


tine  htop 


♦♦♦Output*** 

muzzle  vcl  (m/scc) 

424.10 

max  V  pis  (m/sec) 

4.9 

max  pi  (mpa) 

105.8 

max  p3  (mpa) 

70.4 

max  p5  (mpa) 

87.0 

max  pi  (mpa) 

69.2 

max  pr  (mpa) 

68.6 

max  acc  (k-g) 

3.0 

max  mass  error 

0.03  % 

max  energy  error 

0.02  % 

ballistic  efficiency  = 

31.06  % 

expansion  ratio  = 

10.26 

loss  to  tube  walls  = 

6.82  % 

run  time  =  2.1 

nstep  = 

3429 
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APPENDIX  B: 

INPUT  AND  OUTPUT  FILES  -  SHOT  28 
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Below  is  a  listing  of  the  job  stream  for  a  2  liter  repeatability  shot  28  with  droplet  burning. 
Only  the  parameters  that  have  been  changed  from  Appendix  A  are  described. 

The  injection  and  control  pistons  were  modified  resulting  in  an  initial  liquid  propellant  volume 
which  is  slightly  larger.  The  initial  chamber  volume  is  also  larger  due  to  the  change  in  contour  of 
the  pistons.  The  damper  control  rod  was  redesigned  resulting  in  a  change  in  the  damper  profile. 

A  new  droplet  profile  has  been  derived. 


8.3 


dt28  - 

VIC  -  2  liter  Charge 

10.16 

591.8 

15.5 

43200. 

109000. 

2233.6 

8406. 

641.19 

855.02 

vent4 

17 

-5.4783 

4.8781 

-0.0521 

6.8529 

0.0000 

6.8783 

1.5286 

7.6238 

1.6921 

7.6858 

1.9461 

7.7203 

2.2001 

7.6858 

2.4541 

7.5735 

2.5811 

7.4778 

2.7567 

7.2680 

2.8214 

7.1730 

2.9484 

7.0279 

3.0754 

6.9190 

3.3294 

6.7709 

3.5834 

6.6906 

3.8374 

6.6650 

10.0000 

6.6650 

148.633  0.0 

61400 
no  close 
pisl 


Drop  -  CD=0.95  -  CD5=0.95 

offset  proj  travel 

gun  tube  diameter 

proj  weight 

piston  weight 

vi  v3 

al  a3 

moving  central  bolt 

bolt  radius  versus  piston  travel 


aholc  agres 
rodwt 

pistons  initially  meet  al  point 
piston  resistance 


0.0  0.0 

1.0  0.0 
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disl 

2 


dis.  coeff.  vs.  piston  travel 


dis.  cocff.  vs.  proj  travel  -  tube 


0.0 

487.17 


4 

0.0 

4.0 

3.81 

4.0 

3.82 

0.7 

591.8 

0.7 

5350.0 

9.11 

1.43  5350.0  ' 

4035.5  1.2226 

66.9  .04988 

22.848  .677 

3.45  0.1  0.0 


2 

0.0  .0100 

1000.0  .0100 

prim  3 

450.0  0.004  0.45 

0.0 

heal  1 


steady  state  mass  flux  formulation 
nvo  pth 

iscntropic  flow  into  tube 
proj  resistance 


rhO  kl  k2 
energy  gamma 

surface  tension  kinematic  vi.scosily 
mol  wt  gas  covolumc 
pi  p3  cpsml 


drop4 

droplets 

0.01 

95.2590 

ddr 

pbr 

1.64 

.103 

adri 

bdrl 

1.64 

.103 

adrl 

bdrl 

inject  hot  gas 

primer  ma.ss  injection  lime  heal  loss 
puddle  mass 

no  heat  loss  lo  chamber  walls 


heai2 

300.0  1.0 

shoclc2 

0.1  300. 

1.4  28.84 

tubc4 


bum 

buff2 

139.555  35.04845  0.0 
.95 
912. 

3.45  3.45 

0.8885  1600.  11.3 
0.0  1.4  0.727  0.0 

radius 


15 

0.0000 

3.3147 

0.4724 

3.3147 

0.7264 

3.3122 

0.8534 

3.3084 

2.1488 

3.2652 

2.4028 

3.2603 

2.4714 

3.2634 

2.7254 

3.2680 

2.9794 

3.2736 

3.2334 

3.2802 

3.4874 

3.2878 

3.7414 

3.2964 

3.9954 

3.3061 

4.2494 

3.3170 

4.4527 

3.3264 

2.00C-04 

l.OOc-05 

heat  loss  to  gun  tube  walls 
tube  temp  fudge  factor 
air  shock 
airp  airt 
airgam  airmw 

unsteady  lagrange  distribution 


buffer  -  brayco  783 
apis  ahole  acir 
cd5 
v5 

p5  pout5 

rhO  kl  k2 

eps5  gam5  cv5  b5 

nrat 


line  htop 
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repi 

chaml 


l.OOe-05  l.OOe-08 
10  0 
20.0 


***Output*** 

muzzle  vel  (m/scc) 

411.27 

max  V  pis  (m/scc) 

4.7 

max  pi  (mpa) 

81.9 

max  p3  (mpa) 

56.0 

max  p5  (mpa) 

81.2 

max  pi  (mpa) 

54.6 

max  pr  (mpa) 

54.2 

max  acc  (k-g) 

2.3 

max  mass  error 

0.03  % 

max  energy  error 

0.01  % 

ballistic  efficiency  = 

28.33  % 

expansion  ratio  = 

10.02 

loss  to  tube  walls  = 

7.07  % 

run  time  =  2.3 

nstep  = 

eps  srec 
mf  kwrite 
tmax 

integrate  once 
compute  p3 


3838 
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APPENDIX  C: 

INPUT  AND  OUTPUT  FILES  -  SHOT  51 


S9 


Below  is  a  listing  of  the  job  !>lrcam  for  the  5  liter  Shot  51  with  droplet  burning.  Only  the 
parameters  that  have  been  changed  from  the  Shot  28  model  above  are  described. 

The  initial  liquid  propellant  volume  is  larger  than  the  2  liter  case.  The  volume  is  also  slightly 
smaller  than  that  reported  by  GE  in  Table  2  in  order  to  maintain  consistency  in  the  code  with  the 
given  piston  stroke.  The  damper  vent  area  is  changed.  A  new  droplet  profile  has  been  derived. 
The  projectile  shot  start  pressure  and  resistive  pressure  are  changed  (see  discussion  in  report). 


dt51  -  VIC  -  5  liter  Charge  - 


10.16 

591.8 

15.5 

43200. 

109000. 

5204.86 

8406. 

641.19 

855.02 

veni4 

17 

-5.4783 

4.8781 

-0.0521 

6.8529 

0.0000 

6.8783 

1.5286 

7.6238 

1.6921 

7.6858 

1.9461 

7.7203 

2.2001 

7.6858 

2.4541 

7.5735 

2.5811 

7.4778 

2.7567 

7.2680 

2.8214 

7.1730 

2.9484 

7.0279 

3.0754 

6.9190 

3.3294 

6.7709 

3.5834 

6.6906 

3.8374 

6.6650 

10.0000 

6.6650 

148.633  0.0 

61400 

no  close 

pisi 

2 

0.0 

0.0 

1.0 

0.0 

Drop  •  CD=0.95  -  CD5=0.95 

offset  proj  travel 

gun  tube  diameter 

proj  weight 

piston  weight 

vl  v3 

al  a3 

moving  central  bolt 

bolt  radius  versus  piston  travel 


aholc  agrcs 
rodwt 

pistons  initially  meet  at  point 
piston  resistance 
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disl  dis.  coeff.  vs.  piston  travel 

2 


0.0 

.95 

1.0 

.95 

disl 

dis.  coeff.  vs.  proj  travel  -  tube 

2 

0.0 

1.0 

487.17 

1.0 

fluxl 

steady  slate  mass  flux  formulation 

1 

1.0 

nvo  pth 

flux2 

isentropic  flow  into  tube 

projl 

proj  resistance 

4 

0.0 

10.0 

3.81 

10.0 

4.00 

3.0 

591.8 

3.0 

1.43 

5350.0  9.11 

rhO  kl  k2 

4035.5 

1.2226 

energy  gamma 

66.9 

.04988 

surface  tension  kinematic  viscosity 

22.848 

.677 

mol  wt  gas  covolumc 

3.45 

0.1  0.0 

pi  p3  cpsml 

drop4 

droplets 

0.01 

95.2590 

ddr  pbr 

1.64 

.103 

adrl  bdrl 

1.64 

.103 

adrl  bdrl 

8 

0.0 

.0300 

30.0 

.0300 

40.0 

.0200 

50.0 

.0200 

60.0 

.0100 

70.0 

.0100 

100.0 

.09'’ 5 

1000.0 

.0025 
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prim  3 


inject  hot  gas 

primer  mass  injection  time  heat  loss 


450.0  0.004  0.45 


0.0 

puddle  mass 

heatl 

no  heat  loss  to  chamber  walls 

heai2 

heat  loss  to  gun  tube  walls 

300.0 

1.0 

tube  temp  fudge  factor 

shoclc2 

air  shock 

0.1  300. 

airp  airt 

1.4 

28.84 

airgam  airmw 

tubc4 

unsteady  lagrange  distribution 

bum 

buff2 

buffer  -  brayco  783 

139.555 

35.04845  0.0 

apis  ahole  acir 

.95 

cd5 

1532. 

v5 

3.45 

3.45 

p5  poul5 

0.8885 

1600.  11.3 

rhO  kl  k2 

0.0  1.4 

0.727  0.0 

eps5  gam5  cv5  b5 

radius 

Read  in  damper  rod  radius 

24 

0.0000 

3.3078 

0.0699 

3.3078 

5.1575 

3.1661 

5.6274 

3.1661 

5.7798 

3.1643 

6.0338 

3.1605 

6.2878 

3.1585 

6.3818 

3.1580 

6.6358 

3.1593 

6.8898 

3.1620 

7.1438 

3.1664 

7,3978 

3.1725 

7.6513 

3.1801 

7.9058 

3.1897 
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8.1598 

3.1999 

8.4138 

3.2111 

8.6678 

3.2230 

8.9218 

3.2357 

9.1758 

3.2492 

9.4298 

3.2636 

9.6838 

3.2786 

9.9378 

3.2946 

10.1918 

3.3111 

10.3759 

3.3239 

2.00c -04 

l.OOe-05 

l.OOc-05 

l.OOc-08 

10  0 
20.0 

repl 

chaml 


tine  htop 
cps  srcc 
mf  kwriie 
tmax 

integrate  once 
compute  p3 


***Output*** 


muzzle  vcl  (m/scc) 

612.63 

max  V  pis  (m/scc) 

18.2 

max  pi  (mpa) 

229.1 

max  p3  (mpa) 

149.2 

max  p5  (mpa) 

272.3 

max  pi  (mpa) 

141.2 

max  pr  (mpa) 

139.5 

max  acc  (k-g) 

5.9 

max  mass  error 

0.01  % 

max  energy  error 

0.07  % 

ballistic  efficiency  = 

26.97  % 

expansion  ratio  = 

8.43 

loss  to  tube  walls  = 

5.65  % 

run  time  =  1.9 

nstep  = 

3267 
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APPENDIX  D; 

INPUT  AND  OUTPUT  FILES  -  SHOT  58 


Below  is  a  listing  of  the  job  stream  for  the  5  liter  Shot  58  with  droplet  burning.  Only  the 
parameters  that  have  been  changed  from  the  Shot  51  model  in  Appendix  C  are  described. 

The  major  difference  is  in  the  igniter.  The  igniter  is  a  hybrid  igniter  consisting  of  an  externally 
mounted  solid  propellant  igniter  and  a  liquid  prepositioned  charge  in  the  chamber.  A  new  droplet 
profile  has  been  derived. 
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clt58  -  VIC  -  5  liter  Charge  - 
10.16  591.8 

15.5 
43200. 

109000. 

5204.86  8406. 

641.19  855.02 

vem4 

17 

-5.4783  4.8781 

-0.0521  6.8529 

0.0000  6.8783 

1.5286  7.6238 

1.6921  7.6858 

1.9461  7.7203 

2.2001  7.6858 

2.4541  7.5735 

2.5811  7.4778 

2.7567  7.2680 

2.8214  7.1730 

2.9484  7.0279 

3.0754  6.9190 

3.3294  6.7709 

3.5834  6.6906 

3.8374  6.6650 

10.0000  6.6650 

148.633  0.0 

61400 
no  close 
pisi 
2 

0.0  0.0 

1.0  o.n 

disi 


Drop  -  CD=0.95  -  CD5=0.95 

offset  proj  travel 

gun  tubo  diameter 

proj  weight 

piston  weight 

vl  v3 

al  a3 

moving  central  boll 
bolt  radius  vs.  piston  travel 


aholc  agres 
rodwt 

pi.stons  initially  meet  at  point 
piston  resistance 


dis.  coeff.  vs.  piston  travel 
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0.0  .95 

1.0  .95 

disl  dis.  cocff.  vs.  proj  travel  -  lube 

2 


0.0 

1.0 

487.17 

1.0 

fluxl 

steady  stale  mass  flux  formulation 

1 

1.0 

nvo  pih 

flux2 

iseniropic  flow  into  tube 

projl 

proj  resistance 

4 

0.0  10.0 

3.81  10.0 

4.00  3.0 

591.8  3.0 

1.43  5350.0  9.11  rhO  kl  lc2 

4035.5  1.2226  energy  gamma 

66.9  .04988  surface  tension  kinematic  viscosity 


22.848 

.677 

mol  wt  gas  covolumc 

3.45 

0.1  0.0 

pi  p3  cpsml 

drop4 

droplets 

0.01 

95.2590 

ddr  pbr 

1.64 

.103 

adrl  bdrl 

1.64 

.103 

adrl  bdrl 

8 

0.0 

.0300 

25.0 

.0300 

50.0 

.0100 

75.0 

.0075 

100.0 

.0020 

125.0 

.0010 

150.0 

.0001 

1000.0 

.0001 

prim3  inject  hoi  gas 

280.0  0.004  0.45  primer  ma.ss  injeclion  lime  heal  loss 


10! 


133.0 

hcatl 

hcatZ 

300.0 

shock2 

0.1 

1.4 

tubc4 
bum 
buff2 
139.555 
.95 
1532. 
3.45 
0.8885 
0.0  1.4 

radiu.s 

24 

0.0000 

0.0699 

5.1575 

5.6274 

5.7798 

6.0338 

6.2878 

6.3818 

6.6358 

6.8898 

7.1438 

7.3978 

7.6518 

7.9058 

8.1598 

8,4138 

8.6678 


1.0 

300. 

28.84 


35.04845  0.0 


3.45 

1600.  11.3 
0.727  0.0 


3.3078 
3.3078 
3.1661 
3.1661 
3.1643 
3.1605 
3.1585 
3.1580 
3.1593 
3.1620 
3.1664 
3.1725 
3.1801 
3.1897 
3.1999 
3,21  1  1 
3,22.30 


puddle  mass 

no  heat  loss  to  chamber  walls 
heat  loss  to  gun  tube  walls 
tube  temp  fudge  factor 
air  shock 
airp  airt 
airgam  airmw 

unsteady  lagrange  distribution 

buffer  -  brayco  783 
apis  ahole  acir 
cd5 
v5 

p5  pouts 
rhO  kl  k2 
eps5  gam5  cv5  b5 
Read  in  damper  rod  radius 
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8.9218 

3.2357 

9.1758 

3.2492 

9.4298 

3.2636 

9.6838 

3.2786 

9.9378 

3.2946 

10.1918 

3.3111 

10.3759 

3.3239 

2.00e-04 

1.00c -05 

tine 

htop 

l.OOc-05 

1.00c -08 

ens 

srcc 

10  0 

mf 

kwritc 

20.0 

tmax 

rep  I  integrate  once 

chaml  compute  p3 

***Output*** 


muzzle  vcl  (m/scc) 

617.16 

max  V  pis  (m/scc) 

18.6 

max  pi  (mpa) 

233.1 

max  p3  (mpa) 

151.9 

max  p5  (mpa) 

285.5 

max  pi  (mpa) 

142.8 

max  pr  (mpa) 

143.2 

max  acc  (k-g) 

6.1 

max  mass  error 

0.0!  % 

max  energy  error 

0.07  % 

ballistic  efficiency  = 

27.37  % 

expansion  ratio  = 

8.43 

loss  to  tube  walls  = 

5.70  % 

run  lime  =  2.4 

nstep  = 

4429 
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APPENDIX  E: 

INPUT  AND  OUTPUT  FILES  -  SHOT  65 
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Below  is  a  listing  of  the  job  stream  for  the  7  liter  Shot  65  with  droplet  burning.  Only  the 
parameters  that  have  been  changed  from  the  Shot  58  model  in  Appendix  D  are  described. 

The  initial  liquid  propellant  volume  is  larger.  The  damper  vent  area  is  changed  to  accomodate 
the  longer  piston  stroke.  The  pre-positioned  charge  is  smaller.  A  new  droplet  profile  has  been  derived. 
The  projectile  shot  start  pressure  is  changed  (.see  discussion  in  report). 
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dt65  - 

VIC  -  7  liter  Charge 

10.16 

591.8 

15.5 

43200. 

109000. 

7217.00 

8406. 

641.19 

855.02 

vent4 

17 

-5.4783 

4.8781 

-0.0521 

6.8529 

0.0000 

6.8783 

1.5286 

7.6238 

1.6921 

7.6858 

1.9461 

7.7203 

2.2001 

7.6858 

2.4541 

7.5735 

2.5811 

7.4778 

2.7567 

7.2680 

2.8214 

7.1730 

2.9484 

7.0279 

3.0754 

6.9190 

3.3294 

6.7709 

3.5834 

6.6906 

3.8374 

6.6650 

10.0000 

6.6650 

148.633  0.0 

61400 
no  close 
pisl 
7 


Drop  -  CD=0.95  -  CD5=0.95 

offset  proj  travel 

gun  tube  diameter 

proj  weight 

piston  weight 

vl  v3 

al  a3 

moving  central  bolt 
bolt  radius  vs.  piston  travel 


ahole  agres 
rodwt 

pistons  initially  meet  at  point 
piston  resistance 


0.0  0.0 

1.0  0.0 
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dis.  coeff.  vs.  piston  travel 


disl 
2 

0.0  .95 

1.0  .95 

disl  dis.  coeff.  vs.  proj  travel  -  tube 

2 


0.0 

1.0 

487.17 

I.O 

fluxl 

steady  state  mass  flux  formulation 

1 

1.0 

nvo  pth 

flux2 

iscntropic  flow  into  tube 

projl 

proj  resistance 

4 

0.0 

20.0 

3.81 

20.0 

4.00 

3.0 

591.8 

3.0 

1.43 

5350.0  9.11 

rhO  kl  k2 

4035.5 

1.2226 

energy  gamma 

66.9 

.04988 

surface  tension  kinematic  viscosity 

22.848 

.677 

mol  wt  gas  covolume 

3.45 

0.1  0.0 

pi  p3  epsml 

drop4 

droplets 

0.01 

95.2590 

ddr  pbr 

1.64 

.103 

adrl  bdrl 

1.64 

.103 

adrl  bdrl 

8 

0.0 

.0300 

25.0 

.0300 

50.0 

.0125 

75.0 

.0100 

100.0 

.0060 

125.0 

.0030 

150.0 

.0015 

1000.0 

.0015 
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prim3 

inject  hot  gas 

280.0  0.004  0.45 

primer  mass  injection  time 

115.2 

puddle  mass 

heail 

no  heal  loss  to  chamber  walls 

heai2 

heat  loss  to  gun  tube  walls 

300.0 

1.0 

tube  temp  fudge  factor 

shock2 

air  shock 

0.1  300. 

airp  airt 

1.4  28.84 

airgam  airmw 

tubc4 

unsteady  lagrangc  distribution 

bum 

bufO 

buffer  -  bray  CO  783 

139.555 

35.04845  0.0 

apis  ahole  acir 

.95 

cd5 

1953. 

v5 

3.45 

3.45 

p5  pouL5 

0.8885 

1600.  11.3 

ritO  k  1  k2 

0.0  1.4 

0.727  0.0 

eps5  gam5  cv5  b5 

radius 

Read  in  damper  rod  radius 

23 

0.0000 

3.3103 

5.1753 

3.1661 

9.6203 

3.1661 

9.7727 

3.1643 

10.0267 

3.1605 

10.2807 

3.1585 

10.3746 

3.1580 

10.6286 

3.1593 

10.8826 

3.1620 

11.1366 

3.1664 

1 1 .3906 

3.1725 

11.6446 

3.1801 

n.8986 

3.1897 

12.1526 

3.1999 

no 


12.4066  3.2111 

12.6606  3.2230 

12.9146  3.2357 

13.1686  3.2492 

13.4226  3.2636 

13.6766  3.2786 

13.9306  3.2946 

14.1846  3.3111 
14.387  3.3251 


2.00e-04  l.OOc-05 

l.OOc-05  l.OOc-08 

10  0 
20.0 


rcpl 

ehaml 

***Output*** 

muzzle  vcl  (m/scc) 

707.73 

max  V  pis  (m/scc) 

22.6 

max  pi  (mpa) 

283.7 

max  p3  (mpa) 

184.9 

max  p5  (mpa) 

364.3 

max  pi  (mpa) 

173.8 

max  pr  (mpa) 

169.1 

max  ace  (k-g) 

7.2 

m.ax  mass  error 

0.00  % 

max  energy  error 

0.08  % 

ballistic  efficiency  = 

25.96  % 

expansion  ratio  = 

7.66 

loss  to  tube  walls  = 

5.15  % 

run  time  =  1  5 

nstep  = 

tine  htop 
eps  srcc 
mf  kwrile 
tmax 

integrate  onee 
eompuic  p3 


2555 
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No  of  No  of 

Copies  Organization  Copies  Organization 


1  Office  of  the  Secretary  of  Defense 
OUSD(A) 

Director,  Live  Fire  Testing 
ATTN:  James  F.  O’Bryon 
Washington,  DC  20301-3110 

2  Administrator 

Defense  Technical  Info  Center 
ATTN:  DTIC-DDA 
Cameron  Station 
Alexandria,  VA  22304-6145 

1  HQDA  (SARD-TR) 

WASH  DC  20310-0001 

1  Commander 

US  Army  Materiel  Command 
ATTN:  AMCDRA-ST 
5001  Eisenhower  Avenue 
Alexandria,  VA  22333-0001 

1  Commander 

US  Army  Laboratory  Command 
ATTN:  AMSLC-DL 
Adelphi,  MD  20783-1145 

2  Commander 

US  Army,  ARDEC 
ATTN:  SMCAR-IMI-I 
Picatinny  Arsenal,  NJ  07806-5000 

2  Commander 

US  Army,  ARDEC 
ATTN:  SMCAR-TDC 
Picatinny  Arsenal,  NJ  07806-5000 

1  Director 

Bcnet  Weapons  Laboratory 
US  Army,  ARDEC 
ATTN:  SMCAR-CCB-TL 
Watervliet,  NY  12189-4050 

1  Commander 

US  Army  Armament,  Munitions 
and  Chemical  Command 
ATTN:  SMCAR-ESP-L 
Rock  Island,  IL  61299-5000 

1  Commander 

US  Army  Aviation  Systems  Command 
ATTN:  AMSAV-DACL 
43(X)  Goodfcllow  Blvd. 

St.  Louis,  MO  63120-1798 
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1 


Director 

US  Army  Aviation  Research 
and  Technology  Activity 
ATTN:  SAVRT-R  (Library) 

M/S  219-3 

Ames  Research  Center 
Moffett  Field,  CA  94035-1000 

Commander 

US  Army  Missile  Command 
ATTN;  AMSMI-RD-CS-R  (DOC) 

Redstone  Arsenal,  AL  35898-5010 

Commander 

US  Army  Tank-Automotive  Command 
ATTN:  AMSTA-TSL  (Technical  Library) 
Warren,  Ml  48397-5000 

Director 

US  Army  TRADOC  Analysis  Command 
ATTN:  ATAA-SL 

White  Sands  Missile  Range,  NM  88002-5.s('r 

Commandant 

US  Army  Infantry  School 

ATTN;  ATSH-CD  (Security  Mgr.) 

Fort  Benning,  GA  31905-5660 

Commandant 
US  Army  Infantry  School 
ATTN:  ATSH-CD-CSO-OR 
Fort  Benning,  GA  31905-5660 

Air  Force  Armament  Laboratory 
ATTN:  AFATL/DLODL 
Eglin  AFB,  FL  32542-5000 

Aberdeen  Proving  Ground 

Dir,  USAMSAA 

ATTN;  AMXSY-D 

AMXSY-MP,  H.  Cohen 

Cdr,  USATECOM 

ATTN;  AMSTE-TD 
Cdr,  CRDEC,  AMCCOM 

ATTN:  SMCCR-RSP-A 
SMCCR-MU 
SMCCR-MSl 

Dir,  VLAMO 

ATTN:  AMSLC-VL-D 
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2  Director 

Defense  Advanced  Research  Projects 
Agency 

ATTN:  J.  Lupo 

J.  Richardson 
1400  Wilson  Boulevard 
Arlington,  VA  22209 

4  HQDA 

ATTN:  SARD-ZT,  G.  Singley 
SARD-TT,  I.  Szkrybalo 
SARD-TC,  C.  Church 

B.  Zimmerman 

WASH  DC  20310 

1  HQ,  US  Army  Materiel  Command 

ATTN:  AMCICP-AD,  B.  Dunetz 
5001  Eisenhower  Avenue 
Alexandria,  VA  22333-0001 

1 3  Commander 

US  Army,  ARDEC 

ATTN:  SMCAR-AEE-B,  D.  Downs 

SMCAR-AEE-BR,  B.  Brodman 
W.  Seals 
A.  Bcardell 

SMCAR-AEE-W,  N.  Slagg 
SMCAR-AEE,  A.  Bracuti 
J.  Lannon 
J.  Salo 
D.  Chieu 

SMCAR-FSS-D,  L.  Frauen 
SMCAR-FSA-S,  H.  Liberman 
Picalinny  Arsenal,  NJ  07806-5000 

3  Commander 

US  Army,  ARDEC 
ATTN:  SMCAR-FSS-DA,  Bldg  94 
J.  Fcneck 
R.  Kopman 
J.  Irizarry 

Picatinny  Arsenal,  NJ  07806-5000 

2  Commandant 

US  Army  Field  Artillery  School 
ATTN:  ATSF-CMW 

ATSF-TSM-CN,  J.  Spicer 
Fort  Sill,  OK  73503 
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4  Director 

Benet  Weapons  Laboratory 
US  Army,  ARDEC 
ATTN:  SMCAR-CCB,  L.  Johnson 
SMCAR-CCB-S,  F.  Heiser 
SMCAR-CCB-DS,  E.  Conroy 
A.  Graham 

Watervliet,  NY  12189-4050 

1  Commander 

Materials  Technology  Laboratory 
US  Army  Laboratory  Command 
ATTN:  SLCMT-MCM-SB,  M.  Levy 
Watertown,  MA  02172-0001 

1  Commander,  USACECOM 

R&D  Technical  Library 
ATTN:  ASQNC-ELC-I-T,  Myer  Center 
Fort  Monmouth,  NJ  07703-5301 

1  Commander 

US  Army  Harry  Diamond  Laboratories 
ATTN:  SLCHD-TA-L 
2800  Powder  Mill  Road 
Adelphi,  MD  20783-1145 

1  Commander 

US  Army  Belvoir  Research  and 
Development  Center 
ATTN:  STRBE-WC 
Technical  Library  (Vault)  B-315 
Fort  Belvoir,  VA  22060-5606 

1  Commander 

US  Army  Research  Office 
ATTN:  Technical  Library 
P.O.  Box  12211 

Research  Triangle  Park,  NC  27709-221 1 

1  Commander 

US  Army,  ARDEC 

ATTN:  SMCAR-CCS-C,  T.  Hung 

Picatinny  Arsenal,  NJ  07806-5000 

1  Commandant 

US  Army  Armor  Center 
ATTN:  ATSB-CD-MLD 
Fort  Knox,  KY  40121 
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1 

Commander 

Naval  Surface  Warfare  Center 

ATTN:  D.A.  Wilson,  Code  G31 

Dahlgren,  VA  22448-5000 

2 

Director 

National  Aeronautics  and  Space 
Administration 

ATTN:  MS-603,  Technical  Library 
MS-86,  Dr,  Povinelli 

1 

Commander 

Naval  Surface  Warfare  Center 

ATTN:  Code  G33,  J.  East 

Dahlgren,  VA  22448-5000 

1 

21000  Brookpark  Road 

Lewis  Research  Center 

Cleveland,  OH  44135 

Director 

2 

Commander 

US  Naval  Surface  Warfare  Center 

ATTN:  0.  Dengel 

K.  Thorsted 

Silver  Spring,  MD  20902-5000 

10 

National  Aeronautics  and  Space 
Administration 

Manned  Spacecraft  Center 

Houston,  TX  77058 

Central  Intelligence  Agency 

1 

Commander  (Code  3247) 

Naval  Weapons  Center 

Guns  Systems  Branch 

China  Lake,  CA  93555-6000 

Office  of  Central  Reference 
Dissemination  Branch 

Room  GE-47  HQS 

Washington,  DC  20502 

1 

OSD/SDIO/IST 

ATTN:  Dr.  Len  Caveny 

Pentagon 

Washington,  DC  20301-7100 

1 

Central  Intelligence  Agency 

ATTN:  Joseph  E.  Backofen 

HQ  Room  5F22 

Washington,  DC  20505 

1 

Commandant 

INAFAS 

ATTN:  ATSF-TSM-CN 

Fort  Sill,  OK  73503-5600 

1 

Calspan  Corporation 

ATITJ:  Technical  Library 

P.O.  Box  400 

Buffalo,  NY  14225 

1 

Director 

Jet  Propulsion  Laboratory 

ATTN:  Technical  Library 

4800  Oak  Grove  Drive 

Pasadena,  CA  91109 

8 

General  Elccuic  Ordnance  System 
Division 

ATTN:  J.  Mandzy,  OP43-220 

R.E.  Mayer 

W.  Pasko 

R.  Pate 

1 

IITRI 

ATTN:  Library 

10  W.  35th  Street 

Chicago,  IL  60616 

I.  Magoon 

J.  Scudicrc 

D.  Bair 

L.  Walter 

100  Plastics  Avenue 

1 

Olin  Chemicals  Research 

ATTN:  David  Gavin 

Pittsfield,  MA  01201-3698 

P.O.  Box  586 

Chesire,  CT  06410-0586 

1 

General  Elecuic  Company 
Armament  Systems  Department 
ATTN:  D.  Maher 

2 

Olin  Corporation 

ATTN:  Victor  A.  Corso 

Burlington,  VT  05401 

Dr.  Ronald  L.  Dotson 

24  Science  Park 

New  Haven,  CT  06511 

1 

Honeywell,  Inc. 

ATTN:  R.E.  Tompkins 

MN38-3300 

10400  Yellow  Circle  Drive 
Minnetonka,  MN  55343 
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1  Paul  Gough  Associates,  Inc. 

ATTN:  Dr.  Paul  Gough 
1048  South  Street 
Portsmouth,  NH  03801-5423 

1  Safely  Consulting  Engineer 

ATTN;  Mr.  C.  James  Dahn 
5420  Pearl  Street 
Rosemont,  IL  60018 

1  Sandia  National  Laboratories 

ATTN:  R.  Rychnovsky,  Division  8152 
P.O.  Box  969 

Livermore,  CA  94551-0969 

1  Sandia  National  Laboratories 

ATTN:  S.  Griffiths,  Division  8244 
P.O.  Box  969 

Livermore,  CA  94551-0969 

1  Sandia  National  Laboratories 

ATTN;  R.  Carling,  Division  8357 
P.O.  Box  969 

Livermore,  CA  94551-0969 

1  Science  Applications,  Inc. 

ATTN:  R.  Edelman 
23146  Cumorah  Crest 
Woodland  Hills,  CA  91364 

2  Science  Applications  Inlemaiional 

Corporation 

ATTN;  Dr.  F.T.  Phillips 
Dr.  Fred  Su 

10210  Campus  Point  Drive 
San  Diego,  CA  92121 

I  Science  Applications  Inlemaiional 

Corporation 

ATTN:  Norman  Batiks 
4900  Waters  Edge  Drive 
Suite  255 

Raleigh,  NC  27606 

I  Sundstrand  Aviation  Operations 

ATTN:  Mr.  Owen  Briles 
P.O.  Box  7202 
Rockford,  IL  61125 

I  Vcrilay  Technology,  Inc. 

ATTN;  E.B.  Fisher 
4845  Millersport  Highway 
P.O.  Box  305 

East  Amherst,  NY  14051-0305 
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1  Director 

Applied  Physics  Laboratory 
The  Johns  Hopkins  University 
Johns  Hopkins  Road 
Uurcl,  MD  20707 

2  Director 
CPIA 

The  Johns  Hopkins  University 
ATTN:  T.  Christian 

Technical  Library 
Johns  Hopkins  Road 
Uurel,  MD  20707 

1  University  of  Illinois  at  Chicago 

ATTN;  Professor  Sohail  Murad 
Department  of  Chemical  Engineering 
Box  4348 

Chicago,  IL  60680 

1  University  of  Maryland  at  College  Park 

ATTN:  Professor  Franz  Kasler 
Department  of  Chemistry 
College  Park,  MD  20742 

1  University  of  Missouri  at  Columbia 

ATTN:  Professor  R.  Thompson 
Department  of  Chemistry 
Columbia,  MO  65211 

1  University  of  Michigan 

ATTN:  Professor  Gerard  M.  Faelh 
Department  of  Aerospace  Engineering 
Ann  Arbor,  MI  48109-3796 

1  University  of  Missouri  at  Columbia 

ATTN:  Professor  F.K.  Ross 
Research  Reactor 
Columbia,  MO  65211 

1  University  of  Missouri  at  Kansas  City 

Department  of  Physics 
ATTN:  Professor  R.D.  Murphy 
1 1 10  East  48th  Street 
Kansas  City,  MO  64110-2499 

1  Pennsylvania  State  University 

Department  of  Mechanical  Engineering 
ATTN;  Professor  K.  Kuo 
University  Park,  PA  16802 
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Princeton  Combustion  Research 
Laboratories,  Inc. 

ATTN:  N.A.  Messina 

M.  Summerfield 
4275  U.S.  Highway  One  North 
Monmouth  Junction,  NJ  08852 

University  of  Arkansas 
-Department  of  Chemical  Engineering 
AITTI:  J.  Havens 
227  Engineering  Building 
Fayetteville,  AR  72701 

University  of  Delaware 
Department  of  Chemistry 
ATTN:  Mr.  James  Cronin 

Professor  Thomas  Brill 
Mr.  Peter  Spohn 
Newark,  DE  19711 

University  of  Texas  at  Austin 
Bureau  of  Engineering  Research 
ATTN:  BRC  EME133,  Room  1.100 
H.  Fair 
10100  Burnet  Road 
Austin,  TX  78758 

Jet  Propulsion  Laboratory 
ATTN:  D.  P.  Maynard 

Technical  Group  Supervisor 
Chemical  Processes  Group 
Propulsion  Systems  Section 
4800  Oak  Grove  Drive 
Pasadena,  CA  91109-8099 


117 


No.  of 

Copies  Organization 

1  GS2  Division 

ATTN;  Clive  Woodley 
Building  R31 
RARDE 
Fort  Halstead 

Sevenoaks,  Kent  TNI4  7BT 
England 
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This  Laboratory  undertakes  a  continuing  effort  to  improve  the  quality  of  the  repons  it  publishes. 
Your  comments/answers  to  the  items/questions  below  will  aid  us  in  our  efforts. 

1 .  BRL  Report  Number  __  BRL-TR-315I _ Daig  of  Report  SEPTEMBER  1990 

2.  Date  Report  Received _ 
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